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FROM STOCHASTIC, INDIVIDUAL-BASED MODELS TO THE CANONICAL 
EQUATION OF ADAPTIVE DYNAMICS - IN ONE STEP 

MARTINA BAAR, ANTON BOVIER, AND NICOLAS CHAMPAGNAT 


Abstract. We consider a model for Darwinian evolution in an asexual population with a large 
but non-constant populations size characterized by a natural birth rate, a logistic death rate mod¬ 
elling competition and a probability of mutation at each birth event. In the present paper, we 
study the long-term behavior of the system in the limit of large population (K oo) size, rare 
mutations (u —>■ 0), and small mutational effects (a —■ 0), proving convergence to the canonical 
equation of adaptive dynamics (CEAD). In contrast to earlier works, e.g. by Champagnat and 
Meleard, we take the three limits simultaneously, i.e. u = uk and a — ax, tend to zero with K, 
subject to conditions that ensure that the time-scale of birth and death events remains separated 
from that of successful mutational events. This slows down the dynamics of the microscopic 
system and leads to serious technical difficulties that requires the use of completely different 
methods. In particular, we cannot use the law of large numbers on the diverging time needed for 
fixation to approximate the stochastic system with the corresponding deterministic one. To solve 
this problem we develop a ’’stochastic Euler scheme” based on coupling arguments that allows 
to control the time evolution of the stochastic system over time-scales that diverge with K. 


1. Introduction 

In this paper we study a mieroscopie model for evolution in a population eharacterized by 
a birth rate with a probability of mutation at eaeh event and a logistie death rate, whieh has 
been studied in many works before [|6l|7l|8l|9l[l3l. More precisely, it is a model for an asexual 
population in which each individual’s ability to survive and to reproduce is a function of a one¬ 
dimensional phenotypic trait, such as body size, the age at maturity, or the rate of food intake. 
The evolution acts on the trait distribution and is the consequence of three basic mechanisms: 
heredity, mutation and selection. Heredity passes the traits trough generations, mutation drives 
the variation of the trait values in the population, and selection acts on individuals with different 
traits and is a consequence of competition between the individuals for limited resources or area. 

The model is a generic stochastic individual-based model and belongs to the models of adap¬ 
tive dynamics. In general, adaptive dynamic models aim to study the interplay between ecology 
(viewed as driving selection) and evolution, more precisely, the interplay between the three 
basic mechanisms mentioned above. It tries to develop general tools to study the long time 
evolution of a wide variety of ecological scenarios [|T0l[TTll2T]l . These tools are based on the as¬ 
sumption of separation of ecological and evolutionary time scales and on the notion of invasion 
fitness [fT^IIOl . While the biological theory of adaptive dynamics is based on partly heuristic 
derivations, various aspects of the theory have been derived rigorously over the last years in the 
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context of stochastic, individual-based models [|3|7l[8l|9l[l5l[I3. All of them concern the limit 
when the population size, K, tends to infinity. They either study the separation of ecological 
and evolutionary time scales based on a limit of rare mutations, m —)■ 0, combined with a limit 
of large population miH, the limit of small mutation effects, a —)■ 0, Esin, the stationary 
behavior of the system [fT^ . or the links between individual-based and infinite-population mod¬ 
els [jH. A important concept in the theory of adaptive dynamics is the canonical equation of 
adaptive dynamics (CEAD), introduced by U. Dieckmann and R. Law IfTOl . It is an ODE that 
describes the evolution in time of the expected trait value in a monomorphic population. The 
heuristics leading to the CEAD are based on the biological assumptions of large population and 
rare mutations with small effects and the assumption that no two different traits can coexist. 
(Note that we write sometimes mutation steps instead of effects.) There are mathematically 
rigorous papers that show that the limit of large population combined with rare mutations leads 
to a jump process, the Trait Substitution Sequence, [|6l, and that this jump process converges, in 
the limit of small mutation steps, to the CEAD, @. Since these two limits are applied separately 
and on different time scales, they give no clue about how the biological parameters (population 
size K, probability of mutations u and size of mutation steps a) should compare to ensure that 
the CEAD approximation of the individual-based model is correct. 

The purpose of the present paper is to analyse the situation when the limits of large population 
size, K ^ oo, rare mutations, uk —>■ 0, and small mutation steps, 0, are taken simultane¬ 

ously. We consider populations with monomorphic initial condition, meaning that at time zero 
the population consists only of individuals with the same trait. Then we identify a time-scale 
where evolution can be described as a succession of mutant invasions. To prove convergence to 
the CEAD, we show that if a mutation occurs, the individuals holding this mutant trait can either 
die out or invade the resident population on this time scale, where invasion means that the mu¬ 
tant trait supersedes the resident trait i.e. the individuals with the resident trait become extinct 
after some time. This implies that the population stays essentially monomorphic with a trait 
that evolves in time. We will impose conditions on the mutation rates that imply a separation 
of ecological and evolutionary time scales in the sense that an invading mutant population con¬ 
verges to its ecological equilibrium before a new invading (successful) mutant appears. In order 
to avoid too restrictive hypothesis on the mutation rates, we do, however, allow non-invading 
(unsuccessful) mutation events during this time, in contrast to all earlier works. 

We will see that the combination of the three limits simultaneously, entails some consider¬ 
able technical difficulties. The fact that the mutants have only a iT-dependent small evolution¬ 
ary advantage decelerates the dynamics of the microscopic process such that the time of any 
macroscopic change between resident and mutant diverges with K. This makes it impossible to 
use a law of large numbers as in @ to approximate the stochastic system with the correspond¬ 
ing deterministic system during the time of invasion. Showing that the stochastic system still 
follows in an appropriate sense the corresponding competition Lotka-Volterra system (with in¬ 
dependent coefficients) requires a completely new approach. Developing this approach, which 
can be seen as a rigorous ’’stochastic Euler-scheme”, is the main novelty of the present paper. 
The proof requires methods, based on couplings with discrete time Markov chains combined 
some standard potential theory arguments for the ’’exit from a domain problem” in a moderate 
deviations regime, as well as comparison and convergence results of branching processes. Note 
that since the result of [|6l is already different from classical time scales separations results (cf. 
[fT4l l. our result differs from them a fortiori. Thus, our result can be seen as a rigorous justifi¬ 
cation of the biologically motivated, heuristic assumptions which lead to CEAD. 
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The remainder of this paper is organised as follows. In Seetion and we introduee the 
model and give an overview on previous related results. In Seetion|^we state our results and 
give a detailed outline of the proof. Full details of the proof are presented in the Seetion 
and[^ In the appendix we state and prove several elementary faets that are used throughout the 
proof. 


2. The individual-based model 

In this seetion we introduce the model we analyze. We consider a population of a single asex¬ 
ual species that is composed of a finite number of individuals, each of them characterized by 
a one-dimensional phenotypic trait. The microscopic model is an individual-based model with 
non-linear density-dependence, which has already been studied in ecological or evolutionary 
contexts by many authors lISlI^I^ITSl. 


The trait space X is assumed to be a compact interval of M. We introduce the following 
biological parameters: 

(i) h{x) G M+ is the rate of birth of an individual with trait x E X. 

(ii) d{x) E M+ is the rate of natural death of an individual with trait x E X. 

(iii) iT G N is a parameter which scales the population size. 

(iv) c(x, y)K~^ G M+ is the competition kernel which models the competition pressure felt by 
an individual with trait x E X from an individual with trait y E X. 

(v) UK'm{x) with uk, m{x) E [0,1] is the probability that a mutation occurs at birth from an 
individual with trait x E X, where uk E [0,1] is a scaling parameter. 

(vi) M{x, dh) is the mutation law of the mutational jump h. If the mutant is bom from an 
individual with trait x, then the mutant trait is given hy x + axh G X, where G [0,1] 
is a parameter scaling the size of mutation and his a random variable with law M{x, dh). 
We restrict for simplicity the setting to mutation measures with support included in Z. 

The three scaling parameters of the model are the population size, controlled by the scaling 
parameter K, the mutation probability, controlled by the scaling parameter uk, the mutation 
size, controlled by the scaling parameter ax- The novelty of our approach is that we consider 
the case where all these parameters tend to their limit jointly, more precisely that both uk and 
ax are functions of K and tend to zero as K tends to infinity (subject to certain constraints). 

At any time t we consider a finite number, Nt, of individuals, each of them having a trait 
value Xi{t) G Tf. It is convenient to represent the population state at time t by the rescaled point 
measure, i/^, which depends on K, ux and ax 

1 

i=l 


Let {p, f) denote the integral of a measurable function / with respect to the measure p. Then 
, 1) = NtK~^ and for any x E X, the positive number {of, is called the density of 

trait X at time t. With this notation, an individual with trait x in the population dies due to 
age or competition with rate 



cix,y)iy^{dy). 


d{x) + 


( 2 . 2 ) 
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Let M.[X) denote the set of finite nonnegative measures on X, equipped with the weak topol¬ 
ogy, and define 


M^iX) 



: n > 0, Xi, ...,Xn e 


A’ 


2=1 


(2.3) 


Similar as in ||T3]| . we obtain that the population process, (z//^)i>o, is a A")-valued Markov 
process with infinitesimal generator, defined for any bounded measurable function / from 
M.^{X) to M and for all /i^ G Ai^{X) by 


^ j (j ^ ~ “ UKm{x))b{x) Kii^{dx) (2.4) 

+ j j j - f{fi^)^UKm{x)b{x) M{x,dh) Kii^{dx) 

^ ~ ~ + ^c(a:,r/)/i^(dr/)) K^i^{dx). 

The first and second terms are linear (in jj,^) and describe the births (without and with mutation), 
but the third term is non-linear and describes the deaths due to age or competition. The density- 
dependent non-linearity of the third term models the competition in the population, and hence 
drives the selection process. 


Assumption 1. We will use the following assumptions on the parameters of the model: 

(i) b, d and c are measurable functions, and there exist b, d,c<oo such that 

b{.) <b, d{.) < d and c(. ,.) < c. 

(ii) For all x G A, b{x) — d{x) > 0, and there exists c > 0 such that x G A, c < c(x, x). 

(iii) The support of M(x, .) is a subset of Z n A — x and uniformly bounded for all x G A. 
This means that there exists an A G N such that 

M(x, dh) = Y.t=-APk{x)5k{dh), where Y.t=-APk{x) = 1 for any x G A. 

(iv) 6, d, m G C\X, R) and c G ^^(A^ R) . 


Assumptions (i) and (iii) allow to deduce the existence and uniqueness in law of a process on 
D(M_|_, X4^{X)) with infinitesimal generator (cf. ifT^ ). Note that Assumption (iii) differs 
from the assumptions in [fT3ll because we restrict the setting to mutation measures with support 
included in Z and that it ensures that a mutant trait remains in A. Assumption (ii) prevents the 
population from exploding or becoming extinct too fast. Since A is compact. Assumption (iv) 
ensures that the derivatives of the functions b, c, d and m are uniformly Lipschitz-continuous. 


3. Some notation and previous results 

We start with a theorem, due to N. Fournier and S. Meleard, which describes the behavior of 
the populations process, for fixed u and a, when iT —)■ oo. 

Theorem 3.1 (Theorem 5.3 in ifTSll l. Fix u and a. Let Assumption^hold and assume in addi¬ 
tion that the initial conditions Uq converge for K ^ oo in law and for the weak topology on 
Ad (A) to some deterministic finite measure ^ Ad (A) and that sup^ KfuQ , 1)^] < oo. 

Then for all T > 0, the sequence generated by , converges for K ^ oo in law, in 
D([0, T], Ad(A)), to a deterministic continuous function ^ G ((^([O, T], Ad(A)). This measure¬ 
valued function ^ is the unique solution, satisfying sup^gjo (■Cu 1) < oo. of the integro-differential 
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equation written in its weak form: for all bounded and measurable functions, / : A" —)■ M, 


^t{dx)f{x) = / ^o{dx)f{x)+ ds ^s{dx)um{x)b{x) / M{x,dh)f{x+ah) (3.1) 


'X 


' X 


' X 


! X 


+ ds ^sidx)f{x){{l-um{x))b{x)-d{x)- is{,dy)c{x,y) 


' X 


Without mutation one obtains a convergence to the competitive system of Lotka-Volterra 
equations defined below (see ffT3l ). 


Corollary 3.2 (The special case u = 0 and is n-morphic). If the same assumptions as in 
the theorem above with u = 0 hold and if in addition ,^o = then is given 

by = ^i=i where Zi is the solution of the competitive system of Lotka-Volterra 

equations defined below. 


Definition 3.3. For any (xi, ...,Xn) G we denote by LV{n, (xi, the competitive 

system of Lotka-Volterra equations defined by 


dziif) 

dt 


= Zi i b{xi) - d{xi) - y^c( 

V j=i 



1 < i < n. 


(3.2) 


Next, we introduce the notation of coexisting traits and of invasion fitness (see [|9ll). 


Definition 3.4. We say that the distinct traits x and y coexist if the system LV (2, (x, y)) admits 
an unique non-trivial equilibrium, named z{x, y) E (0, oo)^, which is locally strictly stable in 
the sense that the eigenvalues of the Jacobian matrix of the system LV (2, (x, y)) at z{x, y) are 
all strictly negative. 


The invasion of a single mutant trait in a monomorphic population which is close to its 
equilibrium is governed by its initial growth rate. Therefore, it is convenient to define the 
fitness of a mutant trait by its initial growth rate. 

Definition 3.5. If the resident population has the trait x E X, then we call the following function 
invasion fitness of the mutant trait y 

f{y, x) = b{y) - d{y) - c{y, x)z{x). (3.3) 

Remark 1. The unique strictly stable equilibrium of LV{l,x) is z{x) = hence 

f{x, x) = 0 for all X E X. 

There is a relation between coexistence and invasion fitness (cf. uni). 


Proposition 3.6. There is coexistence in the system LV (2, (x, y)) if and only if 

fix, y) = rix) - c(x, y)z(y) > 0 and fiy, x) = riy) - ciy, x)zix) > 0. (3.4) 

The following convergence result from @ describes the limit behavior of the populations 
process, for fixed a, when K ^ oo and uk 0. More precisely, it says that the rescaled 
individual-based process converges in the sense of finite dimensional distributions to the ’’trait 
substitution sequence” (TSS), if one assumes in addition to Assumption 1 the following ’’Inva¬ 
sion implies fixation” condition. 


Assumption 2. Given any x E X, Lebesgue almost any y E X satisfies one of the following 
conditions: (i) fiy,x) <0 or (ii) fiy,x) > 0 and fix,y) < 0. 
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Note that by Proposition |3.6[ this means that either a mutant eannot invade, or eannot eoexist 
with the resident. 


Theorem 3.7 (Corollary 1 in [[^). Let Assumption^and^hold. Fix a and assume that 

VC > 0, exp(—CiC) <C <C 77 ^—TTTT, as K ^ oo. (3.5) 

K ln(K) 


Fix also X & X and let {N^)k>i be a sequence of'H-valued random variables such that 
{Nq / K) converges for K ^ oo in law to z{x) and is bounded in If for some p > 1. Consider 
the processes generated by with monomorphic initial state {Nq /K)5s^x}- 
Then the sequence of the rescaled processes <^onverges in the sense of finite dimensional 

distributions to the measure-valued process 




(3.6) 


where the X-valued Markov jump process X has initial state Xq = x and infinitesimal genera¬ 
tor 

Afix) = [ (fix + ah) — Mx)) m(x)b(x)z(x) M(x, dh). (3.7) 

Jz 0(x + ah) 


Here we write f{K) g{K) if f{K)/g{K) —)■ 0 when X —)■ cx). Note that, for any 
s < t, the eonvergenee does not hold in law for the Skorokhod topology on D([s, f\,M{X)), 
for any topology Ai(X) sueh that the total mass funetion u i—>■ (z/, 1) is eontinuous, beeause the 
total mass of the limit proeess is a diseontinuous funetion. The main part of the proof of this 
theorem is the study of the invasion of a mutant trait y that has just appeared in a monomorphie 
population with trait x. The invasion ean be divided into three steps. Firstly, as long as the 
mutant population size , !{?;}) is smaller than a fixed small e > 0, the resident population 
size stays close to z{x). Therefore, can be approximated by a linear 

branching process with birth rate b(y) and death rate d{y) + c{y,x)z{x) until it goes extinct 
or reaches e. Secondly, once , l{y}) has reached e, for large K, is close to the solution 
of TC(2, {x,y)) with initial state {z{x),e), w hich reaches the e-neighborhood of (0,V(z/)) in 
finite time. This is a consequence of Corollary]^ Finally, once {of, is close to z{y) and 

{uf, l{a;}) is small, {uf, Ifx}) can be approximated by a subcritical process, which becomes 
extinct a.s. . The time of the first and third step are proportional to ln(X), whereas the time 
of the second step is bounded. Thus, the second inequality in (3.51 guarantees that, with high 
probability, the three steps of invasion are completed before a new mutation occurs. 


Without Assumption]^ it is possible to construct the ’’polymorphic evolution sequence” (PES) 
under additional assumptions on the n-morphic logistic system. This is done in [Q. Finally, 
in Q, the convergence of the TSS with small mutation steps scaled by a to the ’’canonical 
equation of adaptive dynamics” (CEAD) is proved. We indicate the dependence of the TSS of 
the previous Theorem on a with the notation (Xf )t>o. 


Theorem 3.8 (Remark 4.2 in [|9l). If Assumption^is satisfied and the family of initial states of 
the rescaled TSS, Xq , is bounded in If and converges to a random variable X^as a ^ 0, then, 
for each T > 0, the rescaled TSS X)y^2 converges when a ^ oo, in the Skorohod topology on 
D([0, T], X), to the process {xf)t<T with initial state Xq and with deterministic sample paths, 
unique solution of the ordinary differential equation, known as CEAD: 


dxt 

dt 


h [hm{xt) z{xt) dif{xt,Xt)]+M{xt,dh), 


(3.8) 
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where dif denotes the partial derivative of f{x, y) with respect to the first variable x. 


Remark 2. If M(a;, •) is a symmetric measure on Z for all x G X, then the equation (3.81 has 
the classical form, c.f. Ea, 


dxt 1 
dt 2 


m{xt) z{xt) dif{xt,Xt)M{xt,dh), 


(3.9) 


Note that this result does not imply that, applying to the individual-based model first the 
limits {K, uk) (oo, 0) and afterwards the limit a —0 yields its convergence to the CEAD. 
One problem of theses two successive limits is, for example, that the first convergence holds on 
a finite time interval, the second requires to look at the Trait Substitution Sequence on a time 
interval which diverges. Moreover, as already mentioned these two limits give no clue about 
how K, u and a should be compared to ensure that the CEAD approximation is correct. 


4. The main result 


In this section, we present the main result of this paper, namely the convergence to the canon¬ 
ical equation of adaptive dynamics in one step. The time scale on which we control the pop¬ 
ulation process is t/{a\uKK) and corresponds to the combination of the two time scales of 
Theorem |3.7| and 3.8 Since we combine the limits we have to modify the assumptions to obtain 
the convergence. We use in this section the notations and definitions introduced in Section 


Assumption 3. Eor all x G A”, 5i/(x, x) 7^ 0. 

Assumption implies that either Vx G A: dif{x,x) > 0 or Vx G A:9i/(x, x) < 0. 
Therefore coexistence of two traits is not possible. Without loss of generality we can assume 
that, Vx G A, dif{x, x) > 0. In fact, a weaker assumption is sufficient, see Remark |^(iii). 

Theorem 4.1. Assume that Assumptions\Ijand^hold and that there exists a small a > 0 such 
that 


K -C -C 1 and 

exp(-7f") <^uk<^ as 

A In A 


iV —)■ 00 . 


(4.1) 

(4.2) 


Fix Xq G A and let {N^)k>o be a sequence of'H-valued random variables such that NqK~^ 
converges in law, as K ^ 00 , to the positive constant z{xq) and is bounded in L^, for some 

p> 1. 

For each iV > 0, let be the process generated by with monomorphic initial state 
NqK~"^5{xo}- Then, for all T > 0, the sequence of rescaled processes, o<t<T’ 

converges in probability, as K ^ 00 , with respect to the Skorokhod topology on D([0, T], 3W(A)) 
to the measure-valued process 'z{xt)5x^, where {xt)o<t<T A given as a solution of the canonical 
equation of adaptive dynamics. 


dxt 

dt 


h [hm{xt) z{xt) dif{xt,Xt)]+M{xt,dh), 


(4.3) 


with initial condition xq. 

Remark 3. (i) If x* G dX for f > 0, then (4.3) is ^ = 0, i.e. the process stops. 
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(ii) We can prove convergence for a stronger topology. Namely, let us equip the vec¬ 

tor space of signed finite Borel-measures on X, with the following Kantorovich-Rubinstein 
norm: 

ll/iJo = sup { [ fdfit : f G Up-^{x) with sup \f{x)\ < l}, (4.4) 

where Lip^(A’) is the space of Lipschitz continuous functions from A’ to M with Lipschitz 
norm one (cf. ap. 191). Then, for all 5 > 0, we will prove that 


lim P 

if—>oo 


sup 

0<t<T 




> 6 


0 . 


(4.5) 


By Proposition |9 .1 1 this implies convergence in probability with respect to the Skorokhod 
topology. 

(iii) The main result of the paper actually holds under weaker assumptions. More precisely. 
Assumption!^ can be replaced by 

Assumption 3’. The initial state Uq has a.s. (deterministic) support {xq} with xq E X 
satisfying dif{xo, xq) 7 ^ 0 . 


The reason is that since x dif{x,x) is continuous, the Assumption 3 (a) is satisfied 
locally and since x h-)■ dif{x,x) is Lipschitz-continuous, the CEAD never reaches an 
evolutionary singularity (i.e. a value y G X such that dif{y,y) = 0) in finite time. In 
particular, for a fixed T > 0, the CEAD only visits traits in some interval I of X where 
dif{x,x) 7 ^ 0. By modifying the parameters of the model out of / in such a way that 
dif{x, x) 7 ^ 0 everywhere in X, we can apply Thm. |4.1| to this modified process u and 
deduce that yt/KuKry\ has support included in J for f G [0, T] with high probability, and 
hence coincides VtiKuKol. 011 this time interval. 

l + CK 

(iv) The condition uk '^■aK ^hows mutation events during an invasion phase of a mutant 
trait, see below, but ensures that there is no ’’successful” mutational event during this 
phase. 

(v) The fluctuations of the resident population are of order therefore i^-V 2 +“ 

ensures that the sign of the initial growth rate is not influenced by the fluctuations of 
the population size. We will see later that if a mutant trait y appears in a monomorphic 
population with trait x, its initial growth rate is h{y) — d{y) — c{y, x){iyl^, 1) = f{y, x) + 
o{(Tk) = {y- x)dif{x, x) + oiax) since y-x = 0{aK)- 

(vi) exp(iT") is the time the resident population stays with high probability in a 0{eaK)- 
neighborhood of an attractive domain. This can be seen as a moderate derivation result. 
Thus the condition exp(—A'“) uk ensures that the resident population is still in this 
neighborhood when a mutant occurs. 

(vii) The time scale is {KukO'k'^)~^ since the expected time for a mutation event is {Kuk)~^, 

the probability that a mutant invades is of order ax and one needs mutant invasions 

to see a 0(1) change of the resident trait value. This is consistent with the combination of 
of Theorem 13.71 and 13.81 


(viii) Note that the e that we use in the proof of the theorem and in the main idea below will 


not depend on K, but it will converge to zero in the end of the proof of Theorem 4.1 


The constant M introduced below is going to be fixed all the time. It depends only the 
parameters of the model, in particular not on K and e. 
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(ix) The conditions about family of initial states imply that sup^>i sup 4 >o 1)^] < oo 

and therefore, since p > 1, the family of random variable l)}x>i,t>o is uniformly 
integrable (cf. [O Lem. 1). 


4.1. The main idea and the structure of the proof of Theorem 4.1 Under the conditions 


of the theorem the evolution of the population will be described as a succession of mutant 
invasions. We prove that, on the timescale of this result, coexistence of two traits cannot occur, 
namely, when a mutant trait invades the population, the resident trait (i.e the trait that gave birth 
to the mutant trait) dies out. We say the mutant trait fixates in the population. Note that this 
does not prevent coexistence with other mutant traits that do not invade. In order to analyze the 
invasion of a mutant we divide the time until a mutant trait has fixated in the population into 
two phases. 


The two invasion phases: (compare with Figure [T]) First, we prove that, as long as all mutant 
densities are smaller than eax, the resident density stays in an Meux-neighborhood of z{x). 
Note that, because the mutations are rare and the population size is large, the monomorphic 
initial population has time to stabilize in an Meax-neighborhood of this equilibrium z(x) before 
the first mutation occurs. (The time of stabilization is of order and the time where the 

first mutant occurs is of order 1/Kuk)- This allows us to approximate the density of one mutant 
trait yi by a branching process with birth rate h{yi) and death rate d{yi) — c{yi,x)z{x) such 
that we can compute the probability that the density of the mutant trait yi reaches eax, which is 
of order ax, as well as the time it takes to reach this level or to die out. Therefore, the process 
needs 0{a]^) mutation events until there appears a mutant subpopulation which reaches a size 
eax- Such a mutant is called successful mutant and its trait will be the next resident trait. 
(In fact, we can calculate the distribution of the successful mutant trait only on an event with 
probability 1 — e, but we show that on an event of probability 1 — o{aK), this distribution has 
support in {x + axh : h G {1,...,^}}. Therefore, the exact value of the mutant trait is 
unknown with probability e, but the difference of the possible values is only of order ax-) We 
prove in this step also that there are never too many different mutants alive at the same time. 
From all this we deduce that the subpopulation of the successful mutant reaches the density 
eax, before a different successful mutant appears. Note that we cannot use large deviation 
results our time scale as used in dUl to prove this step. Instead, we use some standard potential 
theory and coupling arguments to obtain estimates of moderate deviations needed to prove that 
a successful mutant will appear before the resident density exists a Mecxx-neighborhood of its 
equilibrium. 

Second, we want to prove that if a mutant population with trait yg reaches the size eax, it 
will increase to an Mecri^-neighborhood of its equilibrium density z{ys)- Simultaneously, the 
density of the resident trait decreases to eax and finally dies out. Since the fitness advantage of 
the mutant trait is only of order ax, the dynamics of the population process and the correspond¬ 
ing deterministic system are very slow. (Even if we would start at a macroscopic density e, the 
deterministic system needs a time of order af^ to reach an e-neighborhood of its equilibrium 
density). Thus we can not apply the law of large numbers for density-dependent population pro¬ 
cesses (see Chap. 11 of [fT2l l on our time scales which was used in @ and [0 to approximate 
the population process by the solution of the corresponding competition Lotka-Volterra system. 
This is the main difficulty, which requires entirely new techniques. The method we develop 
to handle this situation can be seen as a rigorous stochastic ’’Euler-Scheme”. Nevertheless, 
the proof contains an idea which is strongly connected with the properties of the deterministic 
dynamical system. Namely, the deterministic system of equations for the case ax = 0 has 
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an invariant manifold of fix points with a vector field independent of ax pointing towards this 
manifold. Turning on a small ax, we therefore expect the stochastic system to stay close to this 
invariant manyfold and to move along it with a speed of order ax- With this method we are able 
prove that, in fact, the mutant density reaches the Mecxx-neighborhood of 2 ( 1 / 5 ) and the resident 
trait dies out. Note that it is possible that a unsuccessful mutant is alive a this time. Therefore, 
we prove that after the resident trait has died out, there is a time where the population consists 
only of one trait, namely the one that had fixed, before the next successful mutant occurs. We 
will divide this phase into several steps. A more detailed outline of the structure of the proof is 
given in Section 


population size 



Figure 1 . Typical evolution of the population during a mutant invasion. 


Note that Figure is only a rough draft and not a ’’real” simulation. 

Notation, (a) If X and Y are two random variables on M, we write X ^ F if we can construct 
a random variable, Y on the probability space as X, such that C{Y) = C{Y), and that for 
all well: X(a;)<F(a;). 

(b) If /i and v are two measures in M.{X), then we write z/ ^ /i if: 

(i) (z/, 1 ) < {/i, 1 ) and 

(ii) sup {x e X : X e Supp (z/)} < inf {x G X : x G Supp (/i)} 

Note that (i) and (ii) imply that, for all / G Lip^(X, [—1,1]) that are monotone increasing 
and for all 0 < f < T, 


f{x)dut < / f{x)dfif 


(4.6) 


IX 


lx 


Convergence: Given T > 0, with the results of the two invasion phases, we will define for all 
e > 0 two measure-valued processes, and in D([0, cxd), A1(X)), such that, for all 

e > 0, 


lim P 

K^oo 


\/t < 


— KuKo'i: 


l,K,t , K , 2,K,t 

l^t ^ l^t ^ 


= 1 , 


(4.7) 
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and, for all e > 0 and i G {1, 2}, 


lim P 

iC—>oo 


sup 

0<t<T/(Kuiicr^) 




= 0 , 


(4.8) 


—)■ 0 when e —)■ 0. This implies (4.5) and therefore the 


for some function 5 such that 5(e) 
theorem. 

More precisely, let 9f be the random time of the ith invasion phase, i.e. the first time after 
9^^ such that a mutant density is larger than eax, and let Rf be the trait of the ith successful 
mutant. Knowing the random variables 9f_i and Rf_i, we are able to approximate 9f^ and Rf: 
After the (i-l)th invasion phase (of the process z/^), we define two random times, 9f'^ and 9f'‘^, 
and two random variables Rf'^ and Rf'^ in X, such that 


lim P 

K^oo 


V^<sup^jGN:0f 


R, 


,K,1 


4 Rf ^ Rf’ and 9f'" ^ 9. 


K 


Thus we define and through 


2 


= forte[9i’ ,9i+i), 


= z^5j^k,2, fort e [0. ’ ,0 .+’i). 



(4.9) 


(4.10) 

(4.11) 


for some appropriate masses z} and zf In fact, zj will be approximately z{Rf'^) for t e 
[9f’^,9f!l) and z^ approximately z{Rf'‘^) for t G [9f’‘^,9ff^). We will prove that the times 
9f'^ and 9f’^ are (approximately) exponentially distributed with parameters of order Kukc^k, 
and that the difference of Rf — Rfi is of order aK- The processes and will 

be constructed by slightly modifying the two processes and in order to make them 
Markovian. This will imply by standard arguments from [|T^ that the processes fiff 2 and 

fiff 2 converge to zixfSx, when ax 0, where Xt is the solution of the canonical equation 
of adaptive dynamics. 

All the remaining sections are devoted to the proof of the Theorem |4.1[ 


5. An augmented process and some eeementary properties 

In the proof of Theorem |4. 1 1 we need to construct a augmented process that keeps track of part 
of the history of the population. In particular, we record the number of mutations that occurred 
before t. 

Let A4^(No x X) = Yfi=i ^9) : ’ 2 , > 0, .^(1),... ,^(u) G No x denote the set of 
finite non-negative point measures on No x A” rescaled by K. We write ^{i) = (^i(i), ^ 2 (*)), 
where ^i(i) G No and ^ 2 (*) G The augmented process, (z>^,L^), is a continuous time 
stochastic process with state space Adp(No x A) x No. The label k of an individual with trait 
(A;, x) denotes that there were k — 1 mutational events in the population before the trait (A:, x) 
appeared for the first time in the population. As in [fT3l . we give a path-wise description of 

Notation. Let ^ Er=i ^ -^f x A) and let = K 

be the number of individuals holding a mutation of label k. Then we rewrite as follows, 

00 an''(At^) 00 

where ^(m''(/x^) 

k=l j=l k=l 


= n. 


(5.1) 
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In fact, the x\^ equal in our situation, because the only variation in the trait 

value is driven by mutational events. We need to define three functions. First, H : A^|;(No x 
X) !-)■ (Mo X A:’)^o is defined as 



/(0,a;?) (0,a;^) . 

(1,E) (1,2:2) • 
( 2 , 0 ;?) (2,o;i) . 

.. (0,o;^o(^)) (0,0) 

(1,2^9111(^4)) 

•• (2, 2;2jj2(^)) (2,0) 

(0,0) 
(1,0) , 
(2,0) , 



i ... . 




Second, h : A1;^(M x 

A) I— )■ (A)^o us given in terms of H by 




hij{^^) = the second component of 


(5.2) 


(5.3) 


i.e., if = {i,x), then hij = x. Third, H : A^^(M x A’) i-)- is defined as follows: 

if Er=i then 


^(/i) ^ (6(a(l)), 6(c^(2)), ..., 6(c^(n)), 0, ...), (5.4) 


where 6(c^(l)) < • • • < 6(c^(«))- 


Definition 5.1. Let (f2,A, P) be a abstract probability space. On this space, we define the 
following independent random elements: 

(i) a A-valued random variable Xq (the random initial trait), 

(ii) a sequence of independent Poisson point measures, (iV^‘^‘**((is, di, d9))k>o, on M+ x M x 
M+ with intensity measure ds En>o ^n{di)dz, 

(iii) a sequence of independent Poisson point measures, {N^'’^^{ds, di, d6))k>o, on M+ x M x 
M+ with intensity measure ds En>o ^n{di)dz, 

(vi) a Poisson point measures, di, d6, dh), on M+ x M x M+ x {—A,... ,A} 

with intensity measure ds En>o ^n{di)dz Eje{-A A} {dh). 

Let = 0 and Uq = then we consider the process defined by the following 

equation 




K T K\ 


(5.5) 






0 J No J M+ 


+ 



^ N'^'’^{ds,di,de) 

^'<c^(^fe,i(^'L))+/NgxA• c(lifc,i(i"A).6)i'^ (<^?)} 

1 


0 JNo ./IR+./{-A,...,A} 




X )+..),!) N-'^^^^-{ds,dt,de,dh). 
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Note that the proeess {uf, Lf)t^o is a Markov proeess with generator 




fc >0 


u + — f{i>, L)j (l — UK'm{x))b{x) Kiy{{k, dx)) 


(5.6) 


+ / (^f {d{x) + J c{x,^2)jy{d0^ Kiy{{k,dx))^ 

+ f [ (f ^(L+i,^aiih) ^ ^ _[_ UKiTi{x)b{x) M{x, dh) Ku{d{k, x)). 

Jnoxa: Jz^ ^ ' ' 

Naturally, the proeess generated by defined in Seetion|^is a projection of the process with 
generator . 

The first elementary property we give is that we there exists a rough upper bound for the total 
mass of the population. 


Lemma 5.2. Under the same assumptions as in Theorem \4.1\ there exists a constant, 1/ > 0, 
such that 

lim P inf{f > 0 : {u^, 1) > 46/c} < exp(V^ii') = 0. (5.7) 

K^oo L J 

Proof. Apply Theorem 2 (a) and then Theorem 3 (c) of @. □ 


6. The First Phase of an Invasion 


In the first phase we show that we can approximate the first time when the density of a 
mutant trait reaches the value eax and the trait value of this mutant trait both on an event with 
probability 1 — o{aK) ■ Such a trait will be called successful. 

Assumption 4. Fix e > 0. Let {R^)k>o be a sequence random variables with values in X. 
Then, there exists a constant M > 0 (independent of e and K) such that for all K large enough 

L^ = 0 and (6.1) 

where e N is a sequence of random variable with \z{R^) — N^kK~^\ < Meax a.s.. We 
call R^ the resident trait. 


Note that Assumption is stronger than the initial condition the assume in Theorem 4.1 
However we obtain with high probability As sumption [Rafter a small time. 

Proposition 6.1. Fix e > 0. Suppose that the assumptions of Theorem |4.i| hold. Then, there 
exists a constant M > 0 (independent from e and K), such that 


lim P 

K^oo 


inf{f > 0 : \{u^, 1) — z{x)\ < Meax} < ln{K)aj/ A inf {f > 0 : L{t) > 1} 


= 1 . 
( 6 . 2 ) 


Since we can assume for the moment that Assumption]^ hold, we do not state the proof here. 
In fact, it can be proven in similar way as Lemma [7^ (a). We begin with several notations, 
which we use in the lemmata below. 


Notation. Fixe > 0. Suppose that Assumption[^[^and|^hold. Let 
and Yj^ G X the trait of the k-th mutant, i.e. 

= inf{f > 0 : Lf = k} and = hk,i 


be the k th mutant time. 



(6.3) 
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We denote by the first time sueh that a mutant density is larger than eax, i-e. 

^kvasion = inf {f > 0 : 3fc e {1,..., Lf} sueh that > eaxK } , (6.4) 


and let i?f be the trait value of the mutant whieh is larger than eaxK at time i.e. 

) with fci = inf |fc > 1 : Dyt^{i>^K ) > eaKK\ . 

^invasion ^invasion J 

Furthermore, let be the first time sueh that [3/a] different traits are present in the 

population, i.e. 


(6.5) 


L^{t) 


^diversity 


inf f > 0 : ^ = [S/a] > , 


( 6 . 6 ) 


k=0 


and similarly let 6*^^ the first time sueh that a ’’2nd generation mutant” oeeurs, i.e. a 
mutant which was born from a mutant born from the resident trait . Note that 


qK 

^mut. of mut. 

Then, we define 


< inf {f > 0 : 3/c G {1,..., Lf} such that \R^ — Y^\ > AaK^) ■ 


= ^i^vasion A ^Lrsity A of mut. A exp(iF° 

The following theorem collects the main results of this section. 


(6.7) 


( 6 . 8 ) 


Theorem 6.2. Fix e > 0. Under the Assumptions^^ there exists a constant M > 0 

(independent of e and K) such that for all K large enough 


(i) z/q = N^kK ^5(0 ,where \z{R^) — N^kK < {M/3)eaK a.s.. 

(ii) We can construct on (17, P) two random variables, Rf'^ and Rf’^, such that 


P 


i?f ^ <Rf < Rf’^ and Rf'^ - Rf'^ < Aa 




-,K,2 


K,2 


,K,1 


K 


P 


_ tdK _ 


= 1 — o(crii:), and (6.9) 
= l-O(e). (6.10) 


The distributions of R^’^ and R^’^ are given in Corollary 


6.10 


(Hi) We can construct on (17, P) two exponential random variables, and with 
parameters of order axUxK, such that 


P 


< E^'^ + ln(iT)cr^^ = 1 - o{aK). 


( 6 . 11 ) 


The distributions ofE^’^ and E^'"^ are given in Lemma 6.7 

Moreover, until the first time of invasion, the resident density stays in an eMax- 

neighborhood ofz(R^), the number of different living mutant traits is bounded by [a/3], and 
there is no mutant of a mutant, with probability 1 — o(ctx). I.e. 


P 


fjK 

^invasion 
= 1 - 


< inf U > 0 : lait°(Ff) 


\Kz{R^)] I > eMaxK] A A 


o{(Tk)- 


( 6 . 12 ) 


Remark 4. The constant M > 0 depends only on a and on the functions b{.),d{.), c(.,.) and 
m(.), but not on K, R^ and e. 


The first lemma in this section concerns the time of exit from an attracting domain. 
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Lemma 6.3. Fix e > 0. Suppose that the assumptions of Theorem \6.2\ hold. Then, there exists 
a constant M > 0 (independent of e and K) such that 


lim aj} ¥ 

K^oo ^ 


inf > 0 : ) - \K^R^)] \ > eMaKK] < 6 


iK 


= 0 . 


(6.13) 


The statement is stronger than the eorresponding one in Thm. 3(e), sinee the diameter 
of the domain eonverges to zero, when K tends to infinity and sinee it eontains a speed of 
eonvergenee to 0 of the probabilities. Therefore, it follows not from the elassieal results about 
the time of exit from an attraetive domain (of. [IT4l l. Our proof is based on a ooupling with a 
disorete Markov ohain and some standard potential theoretioal argument. 


Proof Define 

X* = - \Kz(R^)\ 

and, for all M > 0, 


(6.14) 


Tg = inf{f > 0 : Xt = 0} and = inf{f > 0 : X* > MeaxK}. (6.15) 


Note that Tq and are stopping times with respeet to the natural filtration of X^, whieh is 

equal to a (911°; s < t), and that the prooess (9Jl°)t>o is not markovian. We oan assooiate with 
the oontinuous time prooess Xt a disorete time (non-Markov) prooess Yn whieh reoords the 
sequenoe of values that Xt takes. (This oan be formally defined by introduoing the sequenoes 
Tfc of the stopping times whieh reoord the instanoes when Xt f Xt- and setting = Xt„.) 
Now, we oan oompute 


^[rMeatiK < R ^ 


A nK A nK 

' ' ^diversity ' ' ^mut. of mut. 


(6.16) 


with respeet to the stopping times defined for the disorete time prooess Yn and exploit the natural 
renewal struoture on Y^,. Therefore, we prove the following olaim. 


Claim: For 1 < f X, and K large enough, 

f + fYn = i, Tn-\-l < ^invasion ^ ^diversity ^ ^mut. ofmut.] (6.17) 

< ^ - {c/Ab)K-H + {clb)eaK = p+{i), 

where c, b, c and b are the lower, respectively upper bounds for birth and competition rates. 
Reoall from Remark 1 that the equilibrium z(R^) is equal to ^ and observe that there 

are at most \3/a]eaKK mutant individuals alive at any time t < ^'j^vasion ^ ^^versity ^mut. ofmut.- 
Therefore, for 1 < i C X, and X large enough, 

f + ^\Yn = Tn+1 < ^invasion ^ ^diversity ^ ^mut. ofmut.] (6.18) 

^ {l-ra{R^)uK)b{R’^) 

— (l-m(R^)uK)b{RR+d{R^)+c{R^ ,RK)K-^{\Kz{R‘i)'\+i) 

, d{R^)+c{R^,R^)K-^{\Kz{R^)'\-i)+c\Zla'\eaKK 

^ {\-m{R^)uK)b{RR+d{RK)+c{RTi,RTi)K-^{\Kz{R^)'\-i) 

^ b{R^)-m{R^)uKb{R^) b{R^)-c{R^ ,R^)K-\i-l)+c\3/a']eaKK 

— 2b{R^)-m(RK)uKb{R^)+c{R^,R^)K-^i ^ 2b{R^)-m{RK)uKb{R^)-c(R^,R^)K-^i 

{c/^b)K~^i + (c/6) \3/a]eaK- 

This proves the olaim. Next we introduoe a eoupling, i.e. we define a disorete time proeess Zn 
with the following properties 

(i) Zf) = Yq, 
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(ii) P — Z + 1, P"n+1 — ^ + — -^n — ^n+1 < ^invasion ^diversity ^ ^mut. ofmut.] 

= P [^+1 = i ~\~^yZn = i-, Tn+l < ^invasion ^ ^diversity ^ ^mut. ofmut.] ’ 

(iii) P[z „+1 = Z + 1, Yn+I = i — l\Yn = Zn = i, T„+l < ^invasion ^diversity ^ ^mut. ofmut.] 

= (i) — P[l^n+1 = * + ^\Yn = Z, Tn+l < ^invasion ^ ^diversity ^mut. ofmut.] ’ 

(iv) P \^Zn+l = Z + l|y^ < Zn = Z, Tn+l < ^invasion ^ ^diversity ^mut. ofmut.] ~ P+ (*)’ 

(v) P[Z„+i = i — l\Yn < Zn = Tn+l < ^invasion ^ ^diversity ^mut. ofmut.] = 1 — 


'^mut. of mut. 

. 

invasion 

and the marginal distribution of Zn is a Markov chain with transition probabilities 


P+ 1*1 


Note that by construction a.s. for all n such that A A 9^^^ 


P[Z„+1 = j\Zn = z] = < 


' 1 for z = 0 and j = I 

p+{i) for z > 1 and j = z + 1 

1 — (z) for z > 1 and j = z + 1 

0 else. 


(6.19) 


Now we define a continuous time process, Z, associated to Zn- To do this let (Tj)j^^ be the 
sequence of jump times of Z, i.e. Zt = Zniit e [Tn, Tn+i), defined for all j G N as follows 


T,- - T,_i — 


VP, 


^ — 1 ^ ^invasion ^ ^diversity ^ ^mut. ofmut. 

else, 


( 6 . 20 ) 


where Wj are independent exponential distributed random variables with mean (C'totairate.ff) ^ 


where Ctotairate = 46c(6 + d + c(46c)). By Lemma 
event rate of ,1) and therefore also for 971°. 


5.2 




total rate 


K is an upper bound for the total 


Define = inf{n > 0 : > MeaKK} and Tq = inf{zz > 0 : = 0}. Then, since 


' MeaKK 

Zt > Xt a.s. for all t < 9; 


K 

invasion 


A 0 ^ A 0 ^ 

' ' ^diversity ' ' ^mut. of mut. ’ 


^\_^MeaKK < *"0 ^kvasion ^diversity ^mut. ofmut.] — ^\_'^MeaKK < *'(f] • (6.21) 

yields that, for all M > 32 [S/a] {cb) / {be) such that Zq < ^MeaxK 


9.2 


Applying Proposition 
and large K large enough, 

P [rLaKK < *-0^] < exp . (6.22) 

Next we prove that the process Xt returns many times to zero before it reaches for the first time 
the value MeaKK. More precisely, we obtain first a lower bound for the number of returns by 
the discrete time process Zn- Then we calculate the time for a return to zero. From now on we 
assume that M > 32 [3/a] (c6) /(6c). Define stopping times with respect to the natural filtration 
of Z which records the number of jumps the process Z needs for m zero-returns: 


= inf < zz > 1 : ^ lzi=o = m 


(6.23) 


i=l 


Let = P[r,^ returns < < ^(1 + 1) returns] ^6 the probability that the Markov chain 

returns exactly m times to zero before it reaches the value MeaxK. We have 

Q° = r[Ti„_„^K < r,f] < exp i-M-') , (6.24) 

and, due to the Markov property, for m > 1 , 

Q- = P [ro^ < (1 - Pi [rLaKK < L]r~' Pf [rLaKK < ri] , (6.25) 
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where the last term in the product is smaller than exp Thus, 

Q™ < exp for all m > 0. 


(6.26) 


Let B be the random variable which records the number of zero returns of before reaches 
MeaxK. With other words, B = n if and only if 
obtain that 


Jfl C/WXV/AW 

-Z ^ 


n returns < '^MectkK < '^n + 1 returns WC 


F[B<n] =^Q*< (n + 1) exp . 


(6.27) 


i=0 


Set Ji = T^z and L = T. 


-z — T z for j > 2. For any j, L is the random time 

i it-iuiii j returns (j — l) returns 

between the (j — l)th and the jth zero return of the associated continuous time process Zt and 

B B+1 

Ij < inf{f > 0 : Zt > MeaxK} < (6.28) 

i=l i=l 

We get an upper bound for the probability which we want to compute 


P 


inf <{ f > 0 : ) - \Kz{R^)'\ \ > eM^KK \ < 9 


nK\ 


iK 


< ^P [inf{f > 0 : Zt > MeaKK} < exp(iL") , B = l] +F[B < n]. 


(6.29) 


l=n 


According to (6.28), if i? = I and if in addition more than 1/2 of the I random times Ij in the 
sum are larger than 2l~^ exp(iT"), then inf{f > 0 : Zt > M^eaKK} is larger than exp(iL“). 
Therefore, for all I > n. 


P 


inf{f > 0 : Zt > MeaKK} < exp(iT") , B = I 


(6.30) 


< P 


^ ^ exp(ir°)| ^ ^2 , B I 


i=l 


As mentioned before, Ctotairate.^^ is a upper bound for the total event rate of {uf, 1). Thus we 
can bound the jump times by a sequence of independent, exponential random variables (Vj )jgF^ 
with mean (C'totairate.f^)”^- Namely, 


Tj — Tj_i = Tj — Tj_i )p Vj 


if T < 9^ A 9^ A 9^ 

i — '^invasion '' '^diversity '' '^mut. of mut. ‘ 


(6.31) 


Otherwise the random times T) — T)_i are by definition independent and exponentially dis¬ 
tributed with mean (Ctotairate.?^)”^- The process Z has to make at least two jumps to return to 
zero. Hence, 


^ iTj, for all i eN, 


(6.32) 


where (IFi)igN are independent, exponential random variables with mean (Cjotairate-f^) Thus 


P 


^ l{/j<2i-lexp(ir“)} > z,B — I 


2 = 1 


< p 


El 

2=1 


{W'i<2Z-lexp(ir“)} 


> - 


(6.33) 


Since F\Wi < 21 ^exp(iT“)] = 1 — exp(—Ctotairate-f^^ ^exp(A'")) and {Wi)i>i are indepen¬ 
dent, we obtain that Yl\=i l{Wi<2Z-i exp(A'“)} i® binomially distributed with n = I and p = 1 — 
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exp(—Ctotairate-f^ I ^exp(i^“)). Therefore, the right hand sidc of (6.33) is cqual to 


i 


I . I (l - exp (-Ctotalratei^/ ^ exp(iT“)))* (exp (-Ctotal rate^^/ ^exp(iT“))) *. (6.34) 

i=l/2 

For the following two eomputations we use the elementary faets that (^) < 2* and I < 2\ for 
all / G N and i < 1. We obtain that, for large K enough, the left hand side of ( |6.29[ ) is bounded 
from above by 

( •) (^ “ exp(it:")))* (exp exp(iT")))^"* + F[B <n\ 

l=n i=l /2 


/ 


< ^ - 2' (l - exp (- ^‘°‘Y‘°-^ exp(iT'")))'^' + F[B <n]. 


(6.35) 


l=n 


By (6.271 we see that P [i? < n] = o((Tj^) if the variable n fulfills the following condition 


n -C exp ax- 

Therefore we choose n = [exp {2K°‘)] and get, for large K enough. 


P 


inf <{ f > 0 : ) - \K^R^)] \ > eMaxK \ < 6 


r>K\ 


iK 


(6.36) 

(6.37) 


< ^ 4' (l - exp (-Ctotalrate^/ ^ exp(7f“))) + o(crx) 

/=|'exp(2J^“)] 

oo ^ 

< (^4(1-exp(-Ctotairate^exp(-iF")))'/"j + o{aK) 


Z=[exp(2ir“)] 


< 2 (4^ (1 - exp (-Ctotalrate^ exp(-7f“)))) 


qexp(2X“)l 


+ o(crx) 


<2(42aotalratei^exp(-iT“)) 


irexp(2X“)l 


+ o^ax) 


□ 


< o{Ke ^“) + o{ax), 

where we used that exp(—a:) > 1 — a: for a: > 0 and K exp(iF“") ax- 

In the following lemma we bound Lf, the number of mutants up to time t, from above and 
below by Poisson counting processes. 


Lemma 6.4. Fix e > 0. Suppose that the assumptions of Theorem 6.2 hold and let M be the 
constant of Lemma |6.3| Then, 


lim a 

K^oo 


(^1 _ P V 0 < f < ^ Lf ^ ) = 0, 


(6.38) 


where A^'^ and A^'^ are Poisson counting processes with parameter afuxK and a^uxK 
with 


af = {z{R^) - eMax)b{R^)m{R^), (6.39) 

af = {z{R^) + e{M+ [3/al) ax) (b{R^)m{R^) + Aax) , (6.40) 

and is a constant depending only on the functions 6(.),m(.) and M{.,h) for h G 

{-A •••,41}. 
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Proof. We obtain from the last lemma that 


P 


VO<t<0^ : z{R^) - eMaK < {hA) <z{R^) + e{M+ \2>/a\)aK = 1 - o(cTi^). 

(6.41) 


Therefore, define 




and similarly 



+ J{-A,...,A} 


^{i<K(z{R^)-eMaK), e<b{RK)uKra{R^)M(R^,h)} (6.42) 

X di, d9, dh) 





0 Jfio J {-A,---A} 


'^[i<K{z(RK)+e[M+\l-\yK)} 

^ ^{e<UK {b{R’< )m(R^)M(R^ ,h)+c)’”^’^ A(tk)} 
X N^'^^‘^^^°yds,di,de,dh), 


(6.43) 


Sinee 6^ < 6, 


:)K 


^ i/j^ut of mut ’ mutant trait differs at most A(Tk from the resident trait, R 

we have thatMii'(6(i?'^'''^'' 

tion rate per individual for an appropriate ehoiee of C‘j^ 
proeess with parameter afuKK. By eonstruction, we obtain (6.381. 


K 


Thus, 

m{R^)M{R^, h) + Aax) is a rough upper bound for the muta- 


'tb.rn.M 


. Note that A^’^ are Poisson counting 
' ' □ 


Next we prove that 9Jl^(z/t), the number of offsprings of the kth mutant alive at time t, can be 
approximated by linear birth and death processes. 


Lemma 6.5. Fix e > 0. Suppose that the assumptions of Theorem 6.2 hold and let M be the 
constant in Lemma |6.3| Then, 


lim aj} (1 

K^oo ^ 


P 


^l<k<Lf^,Wt<e^ : zy\t) ^ rnyy) 4 


iK 


K,!/ 


7K,2, 


= 0, (6.44) 


where resp. Z^'‘^{t) are No-valued processes, which are zero until time rjf, the first 

time s.t. 7^ 0, and afterwards linear, continuous time birth and death processes with 

initial state 1 at time and birth rates per individual 

.1 = = b(Y'<) (1 - UKm{\f)) (6,45) 

and death rate per individual 

df’i = diY^) + c{Y^,R^) {z{R^) + MeaK) +c\‘i/a\eaK (6.46) 
resp. df’" = d{Yi^) + c{X^,R^){yR^)-MeaK) . (6.47) 

Furthermore, define Z^'^{t) = Z^'^{Tk + f) and Z^’‘^{t) = Zlf’^{Tk + t), then the processes 
{{Z^’^, Z^'‘^)}k>i are independent and identically distributed. 


Proof. For any t <6^, any individual of 91t*^(r't) gives birth to a new individual with the same 
trait with rate h{Y^) (l - UKmiYy-)) and dies with rate diY^^) + c(Y,y{d^), 
which belongs to the following interval 

d(Yf) + c{Y^^, R^){z{R^)-MeaK), d(Ff) + c(F/, R^){z{R^)+MeaK) + c\3/a]eaK ■ 

(6.48) 
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Thus, let us define, for k < Lqk, 

e<b{Y^){i-UKm(Yk))}^k (ds, di, dO) 





(6.49) 


' Tk JNo >''IR+ 
r'^k+t 



'Tk JNo JM.+ 


{i<Z^’^{s-),e<d{Y^)+c{Yi^,R^){z(Rf()+MeaK)+c\3/a']eaK}^k dO) 


and similarly 


CTk+t 



'Tk JNq JR+ 

rTk+t 



'Tk ./No IR+ 


^{i<Zk'^(s-),d<h(Yk){i-UKm{Yk))]^k {ds,di,d6) (6.50) 

^{i<Zk’^is-), e<d{Yi^)+ciYi^,R‘i)(z{R^)-M€aK)}^k ('^'^5 di, dO), 


and a similar eonstruetion for k > L^k, where the random variables Yj^ are replaeed by i.i.d. 
ones with distribution fx * M{R^, ■), independent of all the previously introdueed random 
variables, where fx is the homothety of ratio ax- Note that, the Poisson point measures 
and are independent of Yjf and and that the proeesses and only depend 

on Y;^ and r^. By eonstruetion, eonditionally on Y,^ = y and = s, the 

proeess Z^’^ is distributed as a linear birth and death proeesses with birthrate h{y){l—ux'm{y)) 
and death rate d{ri) + c{y, R^){z{R^) + Meax) + c\'Z/a\eax, and similarly for Z^’^. In 
particular, the law of {Z^'^ Z^’^) does not depend on r^. Therefore, defining Qk = a{i)t, t < 
Tk, Y^, 1 < ^ < /c — 1), for all bounded measurable functions Fi,... ,Fk on 

D(M+,Z2), 


E 
= E 
= E 
= E 


Fi(Zf’\zf’^)...F,(Zf\zf’^) 

Fi(Zf’\ Zf’^)... Fk-i{Z^^„ Z^^MFk{Z^’\ Zf’^) 
Fi(Zf’\ Zf’^)... Zt\)nFk{Z^^ Z^’^) 

F,{Z^\ Zf’2)... Fk-,{Zt\, ZtD] nFk{Z^\ ZF% 


Qk] 


y^\ 


(6.51) 


where the last equality follows from the fact that the random variable Y^ is independent of 
{Zf'^, Zf’’^) for 1 < £ < /c — 1. Actually, iY^)i<k<L^K random variables, with law 

fx * M{R^, •). This implies by induction that the processes {(Z^’^, Z^’^)}fc>i are i.i.d.. □ 


Let us define Bu = 1 


inf{t>Tfe : rs'}<inf{t>T/f : 9Jt* (!.'t)=0}' 


This random variable indi¬ 


cates if the fc-th mutant population, which appeared at time Tk, invades or not, i.e. reaches 
eaxK individuals before dying out. The following lemma introduces a sequence of i.i.d. ran¬ 
dom variables B^'^) which are 2-tuples of Ber noul li random variables constructed from 


the processes Z^’^it) and Zj^’^{t) defined in Lemma 6.5 such that {B^)k>o is stochastically 


7K,2 


dominated by the sequences (5f'’*)fc> 


fc>0- 


Lemma 6.6. Fix e > 0. Suppose that the assumptions of Theorem \6.2\ hold and let M be the 
constant of Lemma |6.5| Then, 


lim 

K—>-oo 


a 


,-i 

K 


1 -P 


Vl<fc<Lf^ 


Bl'^ ^ B 


K 



0 , 


(6.52) 
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where Bl’^))k>i is a sequence ofi.i.d. 2-tuples of Bernoulli random variables such that 


a.s. Its distribution is characterized by 


aKq^ih) = = 1\ Y,^ = R^ + haK] 


(6.53) 


b{RK) 

0 


Bernoulli 


ifl<h<A, 

otherwise 


and 


axq^ih) = F[Bl’^ = 1\Y,^ = + haK] 


(6.54) 


+ eL'R, 

0 


Bernoulli 


ifl<h<A, 

otherwise, 


where Cl 


and depend only on a, M and Cl (the Lipschitz constant of our param- 


''Bernoulli ' 

eters). Then, for i = 1,2 and k > 1, B]^^ is a Bernoulli random variable of parameter uxpf, 
where 

A 

(6.55) 

h=l 


p1 = 


Remark 5. (i) For all k>l, F[B]2^ = = 1] = 1 — % and is thereby of order e. 


2,K 


Pi 


Pi 


(ii) We use in here the assumption that dif{x, x) > 0 for all x E X. 
Proof Let Z^’^(f) resp. as defined in Lemma 


6.5 


R*’^ — 11 


and define 

for f = 1,2. 


Then, due to the last lemma 


P 


VI < < Lf, : Bl^^ 4 Sf ^ Bl’^ 


= 1 - o(cric). 


For all VI < A; < we obtain with Proposition |9.3| that 


P 


inf{f > Tk : Z^'\t) > eaxK} < mf{f > Tk : Zl^'\f) = 0}| F;, 


KA, 


rK 


\ii,K 

[^fe -R 1h 




(6.56) 

(6.57) 

(6.58) 


= o(exp(-i^“)), 


where, using that f{x, x) = 0 for all x, we have 


bk 


= f(Y,^,R^)-(c(Y,^,R^)M + c[3/a-j)eaK-UKb(Y,^)m(Y,^) (6.59) 

= dJ(R^,R^)(Y,^’^-R^) - (c(Y,^,R^)M + cl3/a])eaK + 0(al), 


and similarly 


b!’^-dt’^=dif(R^,R^)(Y,^-R^) + c(Y,^,R^)MeaK + 0(af) 


j2,K 


yK dKi 


-K 


K\ 


-K r>K\ 


(6.60) 


Reeall that the sequenee {Yjf)k>i used to eonstruet the proeesses Zjf’^ and is a sequenee 
of i.i.d. random variables with distribution M{R^, •). Sinee — ctjf" < 0 if Y^ — R^ < 0, 
we obtain 


pK’ = 1] = E[P[R-|n'"] = 1] 


s E 

/ie{i,...,A} 


dif{R^ ,R^)aKh-(c(Yf ,R^)M+c.\'i/a-\)taK+0{a\) 
b{RK) 


(6.61) 

M{R^,h). 
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Therefore, there exists a constant > 0 (which depends only on a, M and Cl) such that 

(6.62) 


the sum in the right hand side of (6.611 is, term by term, bounded from below by 


and similarly there exists a constant 

= 1] < CTA- 


imciii > 0 such that 


* + cCb...o„„, ) M{R‘\h) 


jK 


(6.63) 


Next, we introduces two couplings, i.e. we define a sequences of i.i.d. 2-tuples of Bernoulli 


random variables ((5^’ , B^’ ))k>i with the following properties 

(i) P [Bl^^ = 0, Bl^ = Q\ Yi^ = rk + Hgk] = P = 0 I Ff = + haK] and 

¥ [Bl^^ = 1, Bl'^ = l\ Yi^ = R^ + haK] = qHh)aK 

(ii) F[Bl’^ = 1, Bl’^ = l\ Y,f = R^ + hax] = F[bI^ = 1 | + hax] 

F[Bl^ = 1 , Bl^ = Q\Yi^ = R^ + haK] = 1 - qHh)cJK. 


and 


By construction, B]^^ < a.s. and B^^ < B^'^ a.s. for all fc > 1 and these random 

variables satisfy ( |6.53[ ) and ( |6.54[ ). □ 


a,K 


R,K 


2,K 


Notation, (a) Fori G {1,2}, define 

rf = inf > 0 : Z^'\Tk + t) = 0 or Z^’^Tk + t)> ecr^ip} . (6.64) 


Obverse that (T^’*) fc>i are i.i.d. random variables that are independent of A^’\ 

(b) Define = inf{/c > 1 : B^ = 1} and = infjfc > 1 : 5^’* = 1}. Then, are 
independent of A^’*, and we have 


F 


r tk,2 


(c) Define R^ = 


K_ 

k 


!}• 


In fact, we prove at the end of this section that F[tjk < 9^] = 1 — o{aK), i-C. Rf is with 
high probability the random variable which gives the value of the next resident trait and tjk , the 
first time where a successful mutant appears, is approximately exponential distributed as stated 
in lemma below. Note that this time is a random time, but not a stopping time. 


Lemma 6.7. Fix e > 0. Suppose that the assumptions of Theorem \6.2\ hold and let M be the 
constant of Lemma |6.3| Then, 


lim aj} 

K^ao ^ 


(^P Tjk < 9 ^ -F [E^’^ 4 Tjk ^ E^’^} H {tjk < 9 ^] ^ = 0, 


( 6 . 66 ) 


where E^’^ and E^'^ are exponential random variables with mean a^pj aKUKK respectively 
0^2 P2 (TkUkK. 


With other words, we have F\E^'‘^ Tjk ^ E^’^\tjk < 9^] = 1 — o((jj^), provided that 
liminf;^_,.oo P[t/k < 9 ^] > 0. 
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Proof. Let be defined as in Lemma 
Then, we obtain by construction, 


and observe that Tjk = inf > 0 : Lf = 



inf {f > 0: 


Tjk ^ inf {f > 0 : Af'^ 


/'"’■}} n |r,« <»''} 


(6.67) 


= P 




o{(Tk)- 


By definition we have that and are geometrically distributed with parameter pf- ax 
resp. P 2 <Jk and A^^^ and ^4^’^ are Poisson counting processes with parameter o^ukK resp. 
u^ukK■ Therefore, the time between each pair of successive events is exponential distributed 
with parameter u^ukK resp. u^ukK. Since the random variables are independent of 
A^’^ and the sum of geometrically distributed many independent, identically exponentially dis¬ 
tributed random variables is exponentially distributed, we get that inf{f > 0 : Af’^ = 
and inf{f > 0 : Af'"^ = are exponentially distributed with parameter a^uxKpf 

respectively a^uxKp^. □ 


In the next lemma we prove that a mutant invades with high probability before the resident 
population exits the neighborhood of this equilibrium, before too many different mutant traits 
are present and before a mutant of a mutant appears. 


Lemma 6.8. Fix e > 0. Suppose that the assumptions of Theorem \6.2\ hold and let M be the 
constant of Lemma \6.3\ Then, 


lim aj/F > ^diversity A exp(iT“) A = 0. 


K^OO 


Proof. We start with proving the following 

P [tsity < a Casion A ofmut.] = o{(Jk)- 

Define 


ryK,,2 / \ _ 

Zk (s) = 


0 


for s < inf{f > 0: Af’"^ = k} 


7K,2 


Tk + s — inf{f > 0: Af’^ = k} ) for s > inf{f > 0:Af’^ = k}. 


K,2 


( 6 . 68 ) 


(6.69) 


(6.70) 


By construction of A^'"^ and the left hand side of (6.69) does not exceed 


P 


k=l 


inf I f > 0 : ^ r3/«l “ ^ < {Kuk(t]+^) 


+ o{aK). (6.71) 


Next, we compute an upper bound for the mutation events that happen before {KukO']^°‘) 
Since A^’‘^ is a Poisson counting process with parameter a^uxK, Chebychev’s inequality im¬ 
plies that 


P 


Ka 


[KukA+‘^) 


_i > 2an a 


K 


K 


< 


Varf 

\ [Kuk( 


.l + c'l 

K ) 


(2 


K —1—a 


r= 


K^-l-a ■ 


at a 


K 


(6.72) 

K 2 

Next we need an upper bound for the lifetimes of the mutants traits, ’ . First, observe that 
the probability that goes extinct after it has reached the value \eaKK~\ converges to zero 
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very fast. More preeisely, Proposition |9.3| and |9.4| (a) imply that 
P infjf > 0 : } < inf|f > = O} < cxd 


(6.73) 


= P 


inf{t>rfc: Z^’^ = 0} < 


CX) 


P 


inf{t>0: [ecTi^iP] } > inf{t>rA:: Z^’^ = 0} 


= o(exp(-7f“)). 


K. 2 

Note that, for eaeh k, , eonditioned on extinetion, is a suberitieal linear birth and death 
proeess (ef. [I18II 1. Let Z^ ’ denote the eonditioned proeess. If Z^ ’ is suberitieal, then eon- 
ditioning has no effeet, otherwise the birth and death rates are exehanged. Denote by 6^ ’ the 


birth rate and the death rate of Z^’^. Then there exist uniform constants, Ci > 0 and 
C 2 > 0, such that Ciax < < C 20 'k, for all k < /^’^. Thus, ^ p. 109 entails, for 

all k < 


7K,2 


P 


Tf< t 






k 

K,2 




o(exp(-it:-“)). 


(6.74) 


The error term o(exp(—if ")) appears since Z^’^, for k < is conditioned on extinc¬ 
tion before reaching the value \eaKK~\ and not only on extinction. Choose t = — 

')-Mn(if), Then, 

dP(l-if) 


P 


^ HK) 


^’^(l-if)-if(dP-6f’^ 

,W,2 

= 1 + 


- o(exp(-if-“)) 


ii\.z £^^,2 

dk -h 


= 1 


b^'\l-K)-Kidr-bt’^) 
0{aKK-^) 


;ir,2 IK, 2 . o(exp(-if “)) 


and hence 


P 


VI < A; < . rj.K ,2 ^ ^CiaK)-^ \n{K)] = 1 - o{aK) 


(6.75) 


(6.76) 


Therefore, we can bound the first summand of (6.71) by times the probability that 

more than [S/a] — 1 mutation events of take place in an interval of length (Cicxi^)”^ In(if). 

More precisely, (|6.71 ) is smaller than 


K -l-a 


202 


P 


^(Cicrx)-! ln{ir) - r3/«l 


+ o((Ti^). 


Thus, for a small enough, the proof of ( |6.69| ) is concluded by the observation that 

\K,2 


P 


A 


(CiaK)-iln(ir) 


> rs/a] 


(6.77) 


(6.78) 


^ g-afnKiV(Ci<7K)-lln(ir) 


E 

i=\Z/a\ — l 


< {a^UKK^CiaK) ^ ln(if)) 


a^UKKiCiax) ^ln(if))* 
rs/ai-i 


l\ 


= Oicr 




where the last equality holds since uxKaj/ ln(if) <C (cri^)“- 
Next, we want to prove that 

P Kut.ofmut. < {Kuk(t]^'")~^ a A Cersity] = o(aic). 


(6.79) 
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Set, for all A > 0, 

G'(A) = E exp J Zq = 1 , (6.80) 

where {Zt,t > 0) is a linear birth and death process with individual birth rate b and indi¬ 
vidual death rate d. Applying the strong Markov property and the branching property at the 
first jump time of Z and using the facts that G(A)^ = E [exp (—A Ztdt) | Zq = 2] and 
E [exp (-Arfc,tj„^p) | Zq = l] = we obtain 

(6.82) 


bG{XY -{b + d + A)G(A) + d = 0. 


Thus, since 


limA^o G(A) = limA^o E [exp (-A /“ Z* dt) 1 ^ = l] 

+ liniAioE [exp (-A/“ Ztdt) | = l] 


= r r. 


extinction 


< oo] + 0, 


which is 0 in the subcritical case and 1 — d/6 in the supercritical case, it follows that 

6 + d + A- 7(6 + d + A)2-46d 


G(A) = 


26 


(6.83) 


Let Z^ ’ (6) = Z^ ’ (xfc -f t), i.e. a linear birth and death process with birth rate 6^ ’ and death 
rate df'’^. Observe that Z^’^(f) dt gives an upper bound for the sum of the lifetimes of all 
individuals with label k. Since the mutation rate of any individual in the population is smaller 
than buK, the probability that an mutant appears, which was born from an unsuccessful mutant 
with label k, is bounded from above by 


K,2 


1 -E 


exp 

< 1 • 


-uxb I zY’‘^(t)dt 


< inf{f > 0 : zY'‘^(t) > eaxK} 


extinction 


(6.84) 


E 


exp ( —uxb 


yK,2 


(t) dt 


T • • <r DO 

' extinction ^ 


+ o(exp(-iL“)). 


7K,2 
‘k 

hand side of 


Since Z^’ (t 


, con ditioned on extinction, is a sub-critcal linear birth and death process, the right 

<^]{uKb) + o(exp(-A'“)) and 




6.84) is equal to 1 — G 

<oc]i^K^ = 




bY^+d.Y^+‘^Kb-Y{bY‘^+dY‘^+UKb)^-‘ibY‘^d 


2b 


K,2 


,K,2 

k 


,K,2 , ,K,2 , T /(jK,2 , rK,2 , tTt . jK,2, 

‘bk +bk -mxfe-V(A -Hfefc’ +UKbY-id^^' b 


K,2 

'k 


2d" 


if 4 >&fc 


• _r iK,2 ^ iK,2 

if^fc >dk 


2bY^+UKb-0(uii<^Y') 

K'" 

2dY^ +UKb-0{uKo]i^) 


4-p ......^ l^?2 

if 4 


• n 7/^,2 . 7 

if^fc >dk 


:K,2 


= 1 - 0{ux(t]Y) = 1 - o{al+^K-^Y- (6.85) 

Note that we used for the second equality that |6fd^’^ | = ^ax for some .^ > 0. By ( 6.72[ ), the 
total number of unsuccessful mutations until is with probability 

1 — o{ax) smaller or equal 2af'(j^^“". Therefore, we finally obtain that the probability to have 
one mutant of a unsuccessful mutant during that time is o{ax)- On the other hand, let Pf" be a 
Poisson counting process with parameter bux^o'xK and (Z/^’^, t > 0) a linear birth and death 


process with initial state 1 and birth rate b ^’i (Y/^) and death rate d^’^{Y /^), then the probability 
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to have one mutant of the sueeessful mutant until the time {KukO']^'^) ^ A A is 

bounded from above by 


P 


PjK.> 7^0 


tURK 


Ao-^ir ^ ^0 


= E 


+ o{aK) 


{P-,, ^o} V'f’-S'pS-p} + *(4S><p) 


( 6 . 86 ) 


_ZK,1 


_ZK.l 


UorK ^ '0 


+ o(aii') 


< (1 — exp{—buK^o'KKtK)) + P 


^K,i 

Ao-jfir 


> fir 




+ o(cri^), 


for eaeh tx, beeause the mutation rate per individual is bounded by bux and there are at most 
eaxK successful mutant individuals alive until 6*i^vasion- we choose tx = then 

by Proposition |9^ all terms in the last line of ( 6.861 are o{ax)- This implies \6.19 ). 

> eaxK'^. Let be a 


Note that we have 6'i^vasion = + inf {f > 0 : 971^ 


IK 


exponential distributed random variable with mean afpf axUxK. Then, 


P 


K 

mut. of mut. 


Tjk + inf {f > 0 : (ur^^+t) > eaxK] > ^^^ersity exp(/f“) A 6, 

> P > {Kuxa],^n~^] - o{ax). 

Let as defined before, then again by Proposition |9.4| 


(6.87) 


P 


Tf/ > \n{K)a- 


-l-a/2 


= P 


eTrK > ln(iT)a^ 


-1-01/2 


^(.gkK < ^0 


= o[(Jx)- ( 6 . 88 ) 


Since ln(iT)(T^^ {KuxO']^°‘) the Markov inequality for the function f{x) = a:", 

where n is smallest even number which is larger than 2/a, yields 

P [E^’^ + Tjr > [Kuxa]^^)-^] < P [E^'^ > {2Kux(j]^^)-^] + o{ax) (6.89) 

/ {2Kux(j]t'"Yn\ 2 ^ 

< , .. .. -—— = 0[ax). 


{a^pfuxKaxY 


□ 


The following lemma shows that there are no two successful mutants during the first phase 
of an invasion. 


Lemma 6.9. Fix e > 0. Suppose that the assumptions of Theorem 6.2 hold and let M be the 
constant of Lemma \6.3\ Then, 


lim (Tf} P 

K-s-oo 


There is a successful mutation in time interval [tjk , 0 


^ . 1 
invasion] 


= 0.(6.90) 


Proof Let mut.(f) the process which recodes the number of successful mutants born after 


r^x until r^x 


P 


t. Then, 

for alH > 0 such that t^k + t < 6^^ : 


K . pK 

succ. mut. 


(f) ^ = 1 - o{ax), 


(6.91) 


where P/^ is Poisson process with parameter a^p^crxUxK. Define Zfjf{t) as in Lemma 
Then P[Vf < 9^ : {i>t) =4 Zfjff)] > 1 — o{ax)- Note that P/^ and Z^’‘^ are independent 


6.5 
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by construction. Therefore, as in the last lemma, or eaeh Ik, 


P 


There is a sueeessful mutation in [tjk , 0{ 


■ 1 

mvasionj 


(6.92) 


< P 


7^ 0 

ecT^f if 


2if,2 

^eaKK < ^0 


+ o(ax) 


< (1 - exp{-a2P2crKUKKtK)) + 




+ o{crK)- 


With tx = by Proposition 9.4, all terms in the last line of (6.92) are oi^ax)- □ 


The following eorollary gives an approximation for the distribution of the next resident trait. 

Corollary 6.10. Fi x e > 0. Suppose that the assumptions of Theorem \6.2\ hold and let M be the 
constant of Lemma 6.3 Then, there exist two X-valued random variables and with 
distribution 


F[RP^ = R^ + axh] = 


K 


M{R^,l)qf{l) 

pf 

M{R^ ,h)qf{h) 


+ 1 — % ifh = l 

P2 


P2 


ifhE {2,..., A}, 


(6.93) 


and 


F[RP^ = R^ + a xh\ = 


M{R^ ,h)qf{h) 

Tif ) 

P2 

M{R^,A)qf{A) 


P2 


+ 1 - J ifh = A, 

C2 


ifh G {1,..., A - 1} 


(6.94) 


such that 


lim ar} 

ir-s-oo 


1 -P 


T)tC,l , r)K j 

"^1 ~v "^1 


^invasion ^ ^diversity ^ ^mut. ofmut. ^ exp(/f 


= 0 . 

(6.95) 


Proof Define 


Rp^ = 


\zK 
rif) 


IK 

7K 


if jir,l _ jK,2^ 


R + ax, otherwise, 


and = 


K,2 _ J ^jK, 

?K 


if/Fl = JK,2^ 


R^ + Aax, otherwise. 


(6.96) 


By eonstruetion of and YjP’\ we have that (6.951 holds. Next, we eompute 


F[YfK ,2 =R^ + axh, = F[Yi^ = R^ + axh, Sf= 1 | Sf= 1] 




jir,2 


F[Yf = R^ + axh, = 1] 

F[bP^ = 1] 

M{R^,h)q^{h) 

P2 


(6.97) 


and P[/^’i ^ IP^] = 1 - = 1 - pf/pf. Since P[i?f ^ = R^ + axh] = 

P[Yi ^.2 = R^ + o'xh, = /^’^] + l{ft,=i}P[J^’^ 7 ^ J^’^] and similarly for Ry'^, we 

deduce ( |6.93[ ) and ( |6.94 ). □ 
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7. The Second Phase oe an Invasion 


Notation. Let us denote 

Cation = inf {t > Cvasion ^ |Supp(i>f)l = 1 ^nd 1) - z{Rf)\ < (M/3)eCTx} (7.1) 

i.e. the first time after Casion ^neh that the population is monomorphie and in the {M/3)eaK- 
neighborhood of the eorresponding equilibrium. 


Again we start with a theorem, whieh summarizes several of the results of this seetion. 


Theorem 7.1. Fix e > 0. Under the As sumptions there exists a constant, M > 0, 

such that, for all K large enough, 

(i) Uq = where — N^kK~^ \ < {M/3)eaK a-s.. 

(ii) At the first time of invasion, the resident density is in an eMaK-neighborhood of 

z{R^), the number of different living mutant traits is bounded by [a/3] and there is no 


mutant of a mutant, with probability 1 — o((Jx). (cf Theorem 6.2 ) 


(Hi) The time between and Cation smaller than 51n(ii')cr^^ with probability 


1 - o{aK)- 

(iv) The trait of the population at time Cation of the mutant whose density was larger 

than eax at time Cva«on’ Supp{u^K ) = , Ri), with probability 1 — o(aK)- The 

fixation 

distribution of Rf can be approximated as in Corollary 
Moreover, until time Cation’ mass of the population stays in the O (a k)- neighborhood 

of z{R^), the number of different living mutant traits is bounded by [a/3], and there is no 
second successful mutant, with probability 1 — o((Tx). 


6.10 


To prove this, we will divide this phase into five steps, illustrated in Figure]^ 

Step 1: From Cvasion ^ Cut. size e’ first time a mutants density reaehes the value e. 
During this period we approximate the mutant density by a eontinuous time branehing 
proeess, whieh is supereritieal (of order ax)- Thus we obtain that ^ is 


of order (In(iF)cr 




Step 2: From ^ to , the first time the mutant density reaehes a value 

(defined in Eq. ( |7.2[ ) below). This step ean be seen as the ’’stoehastie Euler 
seheme”. The idea is that the total mass of the population stays elose to a funetion 
whieh depends only on the density of the sueeessful mutant. This allows to approximate 
the number of mutants by a diserete time Markov ehain until the mutant density has 
inereased by e. Furthermore we eontrol the number of jumps needed to inerease by 
e using upper and lower bounds for one jump time of the assoeiated eontinuous time 
proeess. Then we reeompute the parameters and start again. Iterating, we obtain that 


qK 

'^mut. size 


9, 


K 

mut. size e 


is also of order In (iF)cT 


-1 

K 


Step 3: 


From Cut. sizecj,„ 


until 9 


K 

res. size e 


, the first time sueh that the density of the resident 


trait R^ deereases to the value e. The proof is very similar to the proof of Step 2, the 
only differenee is that we approximate the number of resident individuals by a diserete 
Markov ehain, whieh deereases slowly. 

Step 4: From Cs size e Cs size O’ *^f^® resident trait R^ goes ex- 

tinet. We approximate the dynamies of the resident trait by a eontinuous time branehing 
proeess whieh is suberitieal (of order ax) and therefore goes extinet, a.s., after a time 
of order \n(K)af^. 
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Step 5: From 6 ^^ q until Even if it is unlikely that this time period is larger than 

0, we have to obtain an upper bound for this time. 


population size 



Figure 2. Evolution of the population after the destiny of the sueeessful mu¬ 
tant has reaehed the value ex- 


Notation. Fix e > 0. Suppose that the assumptions of Theorem |7 . 1 1 hold. Set 


and 




oK 

^2 succ. mut. 


inf 

x&X 


-1 


= inf < f > 0 : 1 




> 3 


k=0 


Moreover, for any ^ > 0, 


oK 

^mut. size ^ 


= inf {f > 0 : 3A; > 1 : = \^K]} , 

Csize? = = 


(7.2) 

(7.3) 


(7.4) 

(7.5) 


-1 


and let Sk be a sequenee in K sueh that 1 -C S'x -C ea}^ . 

Remark 6. Using similar arguments as in the proofs of Eemma|6.3[|6.8|and|6.9[ we obtain 


lim (Tt}V 

ir-s-oo 


^invasion “f 5(7 


-1-0/2 ■ 


K 


ln(7f) > ^diversity A ^2 succ. mut. A exp(7f“) 


= 0 . 


(7.6) 


diversity '' 2 succ. mut. A exp(7f“) the total mass of the population 


More preeisely, until the time A 62 ^ 

stays with high probability in the O(aK) neighborhood of z{R^). This ean be proved similarly 
as Eemma |6.3| or |7.2[ Since we have only a approximation of oder ax (not eax), we have 
less precise bounds for the rates of the mutants and for their success probability. Nevertheless, 
we can bound the mutant subpopulations from above by linear branching processes which are 
slightly supercritical of order ax- 
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7.1. Step 1. The following lemma shows that the total mass stays from the beginning (inelud¬ 
ing the first phase) until ^ in the Meax neighborhood of z{x). 


Lemma 7.2. Fix e > 0. Suppose that the assumptions ofTheorem \7.1\ hold. Then, there exists a 
constant M > 0 (independent of e and K) such that 


lim aj} 

K^oo 


inf {f > 0 : |(z/t, 1) — z{R^)\ > Meax] 

< size e A ^2succ. mut. ^ ^diversity ^ exp(7f“) 


(7.7) 


= 0 . 


Proof. The proof of this lemma is very similar to the one of Lemma |6.3[ therefore we omit some 
details. Define 


tdK\ 


Xt = \{hA)K -\Kz{R 


(7.8) 


We assoeiate with the eontinuous time proeess Xt a diserete time (non-Markov) proeess Yn 
whieh reeords the sequenee of values that Xt takes on. 


Claim: For 1 < i < eK, and K large enough, 


— * + l|^n — fTn+l < 


K 


iK 


A 0 

mut. size e ' ' ^2 succ. mut. 


A 1 

' ' ^diversity] 


(7.9) 


< - - [cl^h)K-\ + {2Ci‘^^^AIb)eax ^ p^{i}, 


where is the sum of the Lipschitz constants for the birth, death and competition rate. 

This ean be proven exaetly as in Lemma [6^ using that b{R^) = d{R^)+c{R^, R^)z{R^) and 
that all mutant traits are at a distance of at most 2 Aax from R^, and hence, \b{x) — b{R^) \ < 
C^ax2A, \d{x) — d{R^)\ < Cfax2A and \c{x,y) — c{R^,R^)\ < Clax2A for all traits x 
and y alive in the population. By continuing as in Lemma [O] we obtain ( 1.1) . □ 

tC _ ol let* fV, ^,-11 TZ \ ^ 


Next we prove that ^ 

notation. 


is smaller than ln(7f )ct^^ and we use the following 


Notation. = inf {f > 0 : 1) - z{R^)\ > Meax} A A ^f^ersity 


constant from Lemma \7.2\ Then, 


Lemma 7.3. Fix e > 0. Suppose that the assumptions of Theorem \7.1\ hold. Let M be the 

> {^Lasion + HK)a} 


lim aA ¥ 

K^oo 


9 


K 

mut. size e 


'K 




= 0 . 


(7.10) 


Proof. To prove this lemma we use a coupling with an linear continuous time birth and death 
process. From Phase 1 and the last lemma we know that is with probability 1 — o{ax) 

smaller than 6 ^. Define ki = the label of the first successful mutation. For any t G 
(^imasion) individual of (Fj gives birth to a new individual with the same trait, , 

with rate 


{l-uxmiRf))biR^)e [biR^)-uxb,biR^)], (7.11) 

and dies with rate 

d{Rf)+ [ c{R^,0dM0, (7.12) 

JxxNo 
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which is smaller than dz = li(-Rf) + c{Ri , R^){z{R^) + Meax) + c{e + \3/a]aK)AaK- 

Similarly as in Lemma [O] we eonstruet, by using a standard eoupling argument, a proeesses Zt 
sueh that 


. +t) (7.13) 

invasion ' 

for all t such that ^'j^vasion + ^ ^ 9^ A ini{t > 0 : > eK}. The proeesses Zt is a 

branehing process starting at with birth rate per individual bz = b{Rf-) — bux and 

with death rate per individual dz- For all e < 2 {m+a+i) ’ 

bz - dz > /(7?f, R^) - cax{Me + A{e + \3/a\ax)) > cfx inf (7.14) 

x&X ^ 


Thus Zt is super-eritical of order ax- Let be the first hitting time of level i hy Zt, then by 
Proposition |9.4| 

>T-o^] <exp(-iT“). (7.15) 

Furthermore, we have the following exponential tail bound, see [HI page 41, 


P 


V\eK] < R 

< exp 

(- 



\ L 


emax„<reiri ^n[Tf,x^\^\eK^ < ^o] 


and max„<|-£if] ^n[T^^x'\\'^\eK~\ < ''"(f] — 0{\n{K)ax) (compare with Proposition 
fore. 


(7.16) 
9.3|). There- 


< ln(iF)a^^ > (1 - exp(-a^"^''))(l - exp(-iF")) = 1 - o{ax), (7.17) 
whieh implies the elaim. □ 


—a/‘i\ 


7.2. Step 2. Reeall that the trait of the sueeessful mutant is R^ + a xh where h G {l,...,y4}. 
Due to the regularity assumptions (iv) in Assumption 1, we have the following estimates: 

b{R^ + axh) = b{R^) + b\R^)axh + 0{{axhf) ( 7 . 18 ) 

d{R^ + axh) = d{R^) + d'{R^)axh + 0{{axh)^) ( 7 . 19 ) 

r{R^ + axh) = r{R^) + r'{R^)axh + 0{{axh)'^) ( 7 . 20 ) 

c{R^ + axh, R^) = c{R^, R^) + dic{R^, R^)axh + 0{{axh)^) ( 7 . 21 ) 

c{R^, R^ + axh) = c{R^, R^) + d^ciR^, R^)axh + 0{{axh)^) ( 7 . 22 ) 

c{R^ + axh, R^ + axh) = c{R^, R^) + {dic{R^, R^) + d2c{R^, R^)) axh ( 7 . 23 ) 

+ 0{{axh)^). 


The deterministic system: Although we eannot use a law of large numbers, to understand 
the behavior of the stoehastie system it is useful to look at the properties of the eorresponding 
deterministie Lotka-Volterra system. The limiting system when iT —)■ cx), with ax = 0, takes 
the simple form 


dm^ 

dt 

dm!l^ 


{r{R^) - c{R^, R^){m^t + , 


(7.24) 


dt 


{r{R^) - c{R^,R^){w!l + m^l^)) . 


(7.25) 
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Figure 3. Vector field of the unperturbed system 


The corresponding vector field is depicted in Figure This system has an invariant manifold 
made of fixed points given by the roots of the equation 

= r{R^)/c{R^, R^) = z{R^), (7.26) 

with > 0. This manifold connects the fixed points of the monomorphic equations, 

{z{R^),0) and {0,z{R^)). Note that z{R^) has the interpretation of the total mass of the 
population in equilibrium. A simple computation shows that the Hessian matrix on the invariant 
manifold is given by 




(7.27) 


The corresponding eigenvectors are (1,-1) with eigenvalue 0, and (m°,f(i?^) — mP) with 
eigenvalue —c{R^, R^)z{R^). 

It follows that the perturbed system 

{r{R^) - c{R^, R^)rrPl - c{R^, R^ + aKh)nPp) (7.28) 


dt 

dt 


= m!l^ {r{R^Raxh) — c{R^+aKh^R^)rrPl — c{R^Raxh^R^ +aKh)rrLl'^^^ 


has an invariant manifold connecting its fix points {z{R^), 0) and (0, z{R^ + axh)), where 
z{R^ + axh) = r{R^ + aKh)/c{R^ + axh^R^ + aKh) in an ciK-neighborhood of the 
unperturbed invariant manifold (see Figure]^. Thus the perturbed deterministic system will 
move quickly towards a small neighborhood of this invariant manifold and then move slowly 
with speed 0{aK) along it. Since the invariant manifold is close to the curve rrp + rrP'^ = 
z{R^), it is reasonable to choose as variables Mt = The motion of the system will 

then be close to the curve ) defined by the condition that the derivative of Mt vanishes for 
Mt = Since 

^ = Mt{r{R'^)-c{R^,R^)Mt) (7.29) 

- [{dic{R^, R^) + d 2 c{R^, R^)) Mt - r'{R^)\ aKhm^^ + 0(a|). 

Setting the right hand side to zero yields the leading orders in ax 

V'(i?^) {dic{R^, R^) + d2c{R^, R^)) 


r{R 


K\ 


c{R\R’^) 


p{mt^) = z{R^) + aKhm!p 


+ 0(4). (7.30) 
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We expect that the stochastic system also evolves along this curve. I.e., we will show that 


increases while the total mass stays close to the curve defined in (7.30). 
Define the function 


(t){y) = z{R^) + aKhy 


r\R^) dic{R^, R^) + d 2 c{R^, R^) 


r{R^) 


c{R\R^) 


and the stopping time 

^Lr^q) = - Ct.sizei(e/2) ^ \{^t, 1) “ 0(^(e/2))| < {M/3)€aK} 


(7.31) 


(7.32) 


The dependence of 0 with respect to the mutant density allows us to decompose the increase of 
the mutant density into successive steps during which the total mass does not move more than 

Meax- 

Lemma 7.4. Fix e > 0. Suppose that the assumptions of Theorem ]/. l\ hold. Then, there exists 
a constant M > 0 (independent ofe, K and i) such that and for all 2<i<2e 

(a) Soon after sizei{ei 2 y population size is close to 0(f |).' 


lim a f}¥ 

K^ao ^ 


^near(f>(i^) ^ (,^mut. sizei{el2) F ^2 race. mut. ^ ^diversity 


(7.33) 


= 0 . 


A inf {t > C. .ize^ie/ 2 ) : 3/c > 1 : fm^(z>,) = [(^ ± l){e/2)K]} 

(b) A change of order efor the mutant density takes more than o(^af^) time: 

lim P[inf {t > C. .ize^i./2) : 3^ > 1 : m\h) = [(* ± l){e/2)K] } (7.34) 


K^oo 


^ sizei{e/2) ^ ^near(j)(i^) ^ ^2succ. mut. ^ ^diversity 


= 0 . 


(c) At the time when the mutant density has changed of order e the total population size is still 
close to 


lim ad ^ 

ir-s-oo ^ 


inf -0(i(e/2))| > Meax} <9. 


^ k9^ (1 331 

2 succ. mut .' ' ^diversity ) 


A inf {( > C. ,,„(,/2) : 3*: > 1 : = [(i ± l)(£/2)ii'l } 


= 0 . 
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(d) A change of order efor the mutant density takes no more than (iax) 


- 1 - 0/2 


time: 


lim (Jj} ¥ 


GLt.size{i+l){e/2) > ^ A 6. 


K 

2 succ. mut. 


A Osversity A inf ^ ^ 1) “ 0(^(e/2))| > McaK 


iK 


(7.36) 
= 0 . 


Remark 1. For each e > 0, Lemma [T!4| implies that the mutant density reaches the value 
with high probability, since e is independent of K. Moreover, for all e > 0, 


P 


air 


>{6 

and 


K 

mut. size e 


+ ln(iF)a"^"“/') A 0, 


air 


A9fz 


K 


= o^ax) 


P 


I 


> 1) - <^(C'cross)l > Meax = o{ax). 


(7.37) 

(7.38) 


Proof. We will prove the lemma by induction over i. Base clause: Compare with Lemma 7.2 
and |7. 3 1 that there exists a constant M > 0 such that I iunK , 1) — 0(0) I is smaller than Meax 
and that A both with probability 1 - o{(Jx)- 

Induction step form f — 1 to f: Assume that the lemma holds true for i — 1, then be prove 
separately that (a)-(d) are true for i, as long as f 

Proof of (a) for i by assuming that the lemma holds for i — 1. In the proof we use the following 
notation 


qK ^ qK 
^i — ^2 succ. mut. 


A A mf{« > ,, : 3A^ > 1 : OT*(S,) = ^(i ± l)(€/2)A'l}. 

(7.39) 

Note that 9^ differs from 9^ defined in Lemma 7.3 We will prove (a) provided it happens 


before 9f and we use the estimates of step (b) for i to prove that it indeed happens before 0/ 
with high probability. 

If the Lemma is true for z — 1, we know that (with (d)) 


K 


P 


■( /2i’ ~ l)(e/2))l < Meax 

mut. size z(€/2) 


= 1 - o{ax)- 


Since (j){x) — (j){y) = 0{h{x — y)ax), we have with probability 1 — o{ax) either 
inf {t > Cut. size 7 ( 6 / 2 ) : - 0(z(e/2))| < {M/3)eax} = 9, 

which implies (a) for z, or at least 


K 

mut. size i{^/2) ’ 


I ■ ■( /2) ’ “ '^(*(^/ 2 )) \<[M + 

mut. size z(e/2) ' 


h 


r'{R^) 

r{R^) 


idic{R^ ,R^)+d2c{R^ ,R^)) 
c(R'=,iZfc) 


(7.40) 

(7.41) 

(7.42) 


Similarly as in many previous lemmata we want to couple K{ut,l) with a discrete time Markov 
Chain. Therefore, let 

Xl = \K{i>„l)-\ff{e/2))K]\, (7.43) 


andT* = Cut. size 7 ( 6 / 2 ) 


jK 


and (T^)fc>i be the sequences of the jump times of (z/^, 1) after Cut size i{e/ 2 )- 
Then let be the associated discrete time process which records the values that XI takes after 

time Cut. size 7(6/2)- 


Claim: There exists a constant C^fekhative 


> 0 such that for all \C^dtrtative^^KK~\ < j < \eK] 


and K large enough, 


PK+1 = 3 + 1|A = 7. r„+i < A'] < 2 - =: p/. 


(7.44) 


Moreover, we can choose 

CdfrLtive = sup ^Wx) + “ dic{x,x) - d 2 c{x,x)\). (7.45) 




r(x) 
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If (z/j, l)i^ > \4){i{e/2))K~\ at time t = T^, then 1)K = \(f>{i{e/2))K] + and eon- 
ditionally on the left hand side of (7.44) is equal to the probability that the next event is a 
birth. Namely, 




Sfc>0 + d{hk,i{yTi^)) + c(hk,i(i>Ti), Odi^TiiO j 

< [Hdi'^)Ek>o^'(^n) + aKhb'(RE^^H^n) + Cl2AaK[3/a-laK€K + 0(a^j,K)] 
X E.>o (Mfl'") + £((«*■) + E»0 


(7.46) 




+ axWi(i>TA) & ) + a(^ ) + 


K 


(CE’")2AaK \3la\aKeK - 0{a\K) 


-1 


For the inequality we have used the faet that, eonditioned on T„ < Qf, there at most aK^\3/a\ 
many unsuccessful mutant individuals which differ at most 2AaK from the resident trait . 
Since Ek>o^^i^T{) = (^r^, 1)-^ which equals [0(i(e/2))iF’| + j conditioned on j = Y^, the 
right hand side of the last inequality is smaller or equals 


h{RE + aKhh\RE + 0(4) 


X 


h{RE + d{RE + c(i?^, RE 


r<A(i(e/2))ir]+i IJ + a j + 


ir 


and by definition of 0 the denominator equals 

jK 


gjffci (r, , \ 

2 b{REaK + R^i - 0{al) 


axh 


*(^/2) ( + d2c{R^,R^) ] + 


2 


0 ( 4 ) 

(7.47) 


(7.48) 


X 


- b'iRE + d\RE + 


Thus, we obtain that the right hand side of (7.46) is bounded from above by 

t{e/2) (^ - d,c{R^, RE - d 2 c{R^, R^ 


1 _ c(Rl^,Rl^) -j.-I _ a^h 

2 3b(r?^) 


4b(iJ-K) 


(7.49) 


+ 


r</)(i(e/2)).ZV] +j 




+ 0 ( 4 ). 


In the case where {vt, l)7f < [0(i(e/2))iT] at time t = T^, we obtain the same inequality but 
with an opposite sign in front of the third term. Since 


r'{R^) an'^+7/yi) r'{R^)K 


2 +i?^) 


K +((e 


^-;^-{d,c{RERE + d2c{RERE) 


SOT'^i (u„i) 

_ n 

K 


<(e/2) T^-d,c{R^,R^)-d 2 c{R^,R^) 


K t3K\ 


?K DX^ 


(7.50) 
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we deduce the claim. Since we choose M such that M > we can construct a Markov 

chain such that > Y^, a.s., for all n such that < Of A inf{t > size i(e/ 2 ) • 
1) — 0(i(e/2))| < ^MeaK} and the marginal distribution of Z„ is a Markov chain with 

Zn = Yn and transition probabilities ^ „ 

° ° (pf forji>landj 2 =Ji + l, 

P[^n+i = j 2 \Zn = Ji] = ll-pf for ji > 1 and ja = ji + 1, (7.51) 

0 else. 


Let C*exit — 2z4 


r'(x) {dic{x,x)+d2c{x,x)) 

r{x) c(x,x) 


. Then,by applying Proposition 


obtain, for all a < (M + C'exit)eai^iT and large K large enough. 




9.5 


(b), we 
7.52) 


inf{n > 0 : Z; > 2(M + Cexit)ecri^iT} < inf{n > 0 : Z;, < {M/3)eaKK} 

< exp {—K°‘). 

Next define S* = infjn > 0 : Z* < ^MeaxK}. This is the random variable, which counts 
the number of jumps Z* makes until it is smaller than eaxK. Note that — Tf), the times 
between two jumps of Xf are exponential distributed with a parameter {h{R^) + d{R^) + 
c{R^, R^)z{R^))z{R^)K + 0{aKK), if is smaller than Of. Thus, 

(7)Vi - Ti) 4 El (7.53) 

where {El)i>o are i.i.d. exponential random variables with parameter infajg;^ b{x)z{x)K. There¬ 
fore, 

^neariji(i|) ^ ^mut. sizei(e/2) ^ 0^ 


P 


< P 


"2 / 


(7.54) 




1=0 


+ p 


^ ^mut. sizei(e/2) ^ ^near0(i|) 


K 


Our next goal is to find a number, n*, such that F[B^ > rii] is o(aj^). Since the transition 
probabilities of Z* do not depend on the present state, we have that Z* — Zq has the same law 
as ELi where (I 4 O fcgN is a sequence of i.i.d. random variables with 

F[V^ = l]=pf and P[f^,* = -l] = l-p^ 


and E 


= —2eaK and |V)!| = 1. Furthermore, we get 
P [B' <ni] > P [inf {j > 0 : Zj - Zq < -[(|M + Cexit)ecr^7T] } < n*] 


(7.55) 

(7.56) 


> P 


rii 




.k=l 


and by applying the 

Hoeffding’s Inequality: (Appendix 2 in Il22l ): Let Fi,..., be independent random vari¬ 
ables, X > 0 and Oj <Yj — E[y)] < bj are bounded for all j. Then 

(7.57) 


P 


® K'] > a: < exp [-2x^ ~ 

L j = l -I j = l 


we obatin 


P 


Y E* > -2eaKni + (rii) 
.k=l 


\l/2+a/2 


< 2exp(-(ni)“). 


(7.58) 
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With Hi = \K{\M + Cexit)], we get -2eaKni + < -[(fM + 

A'“ 2 +“ ^ aK- Applying the exponential Chebyehev inequality (with A = A'") 


sinee 





/ r.K'dM+Cex.,)! \ ■ 

p 


< exp(—AS'i^)E 

exp A ^ El 


i=0 


[ \ '=0 /J 


(7.59) 


inf^^x b{x)z{x)K 


rx(§M+c,,„)i+i 


<exp(-A^i^), u \-( \ 

ml^(zxb[x)z[x)K - A 

< exp (^—XSk + (rA"(|Af + Cexit)] + 1 ) lu + 

|M + Cexit + 1 


inf^^x b(x)z(x)K-\ 


< exp — XSk + A 


ini^^xb{x)z{x) 


+ 0{X^K-^) < exp (-i^“) 


Henee, the left hand side of (7.54) is bounded from above by 

exp(-A-“)+2exp(-(/f(|Af+C„i,)r)+n* [«f < («L. + Sx) A 

This proves the lemma, if we ean show that 

< (Ct. sizei(./2) + Sk) a 0^ar0(if)] = 

Aeeording to Remarkand Lemma |73| we have that 

^ [^2 succ. mut. A ^diversity A ^nput. sizei(e/2) o((J^). 

Therefore, the following proof of (b) for i implies (a) for i. 


(7.60) 

(7.61) 

(7.62) 

□ 


Proof, of (b) for i by assuming that the lemma holds for i — 1. Note that the random elements 
S®, T\ L®, IL®, X®, F® and Z® are not the ones of the last proof. They will be defined during 
this proof. In faet, the strueture of the proof is similar to the one of (a), exeept that we prove a 
lower bound for the time of a ehange of oder e for the mutant density instead of upper bound for 
the time of a ehange of oder eaK of the total mass. We eouple IDTf®, for t > 9^^^ sizei(e/ 2 )’ ^ 
diserete time Markov ehain (depending on i). Therefore, let Tq® = 9^^^ size j(e/ 2 ) i^k)k>i 
the sequenees of jump times of OJtf® after 9^^^^ sizei(e/ 2 )- Furthermore, let (F^®)„>o be the diserete 
time proeess whieh reeords the values that takes i.e. Yq = 9Il^^(z/jn) = \Ki{e/2)] and 
{brp^). Observe that if 


nK . nK 


I) Ainf{f > 0n,ut.sizei(e/2) ^ \{h,t) - (j){i{e/2))\ > 2(M + Cexit)eai^A:}, (7.63) 


we know from the inequality (7.52) that the probability that larger than inf{t > 

mut. sizei(e/ 2 ) ^ 1) “ 0(f(e/2))| > 2{M + Cexit) 6 cr^-iT} is Smaller than exp(-A'“). Define 


9 ^ 




inf{f > 9^^, size*(e/2) ^ K^t, 1) “ 0(*(e/2))| > 2(M + Cexitjea^AT} (7.64) 

A A R^ A R^ 

' ^near(^(i|) ' ' ^2 succ. mut. ' ' ^diversity 
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and Cfitness = inf^eA’ dif{x, x )/ b. Then, for all — [ fK~\ < j < [ ^K] , for K large enough and 
for e small enough, we have that 


P 


r:+. = +] + i|i;‘ = 


K 


(7.65) 


2 P 2^fi*^ness 2 P 274Cyfltness 


sinee the left hand side of (7.65) is equal to the expeetation of the probability that the next event 
is a birth without mutation eonditioned on . Namely, 


h{R^ + aKh){l — UKm{R^ — axh)) 


b{R^ + aKh) + d{R^ + axh) + fj>^^;^;c(R^ + crKh, ^)diyT„(0 


(7.66) 


= b{R^ + axh) b{R^ + axh) + d{R^ + axh) + c{R^ + axh, R^) ( 0(i|) 


■)K 


K 


jK 






K 


+ c{R^ axh, R^ + hax) T2(Af + Cexit))) P 0{uk) 

= b{R^+ axh) 2b{R^+ axh) - f{R^+ axh, R^) + c{R^ + axh, R^) (^0(i|) - /rk^rk) '^ 
+ aMciR^, R^) P 6 {ecTKCl{\ f] P 2(M + Cexit)))] P 0{uk). 


for some G (—1,1). By definition of 0 of (7.66) is equal to 

b{R^ + aKh) \2b{R^ + aRh) - dif{R^, R^)aKh + c{R^ + axh, R^)<TKh (7.67) 

^ +^i {^(^kCI ([f] + 2(M + Cexit))) 

P P cr^ + Ur) 

= b{R^ + aRh) \^b{R^ + OKh) - aRh (l - d,f{R^, R^) 


-1 


P ^i(ecr_fs'C'£( [-] + 2(M + Cexit))) P P P '^k) 


1 I ■ec{R’^,R’^ 


ec{R^,R^) \ dif{R^,R^) , . (C|Xr|]+2(M+Cexit)) , , _2 


b{RK) 


P eaR^i- 


b{RK) 


P P P '^k)- 


■ e c{R^,R^) 


'r(R'P) ' ^ obtain (|7.65l). Thus we 


Then, beeause i < 2e ^Ccross implies that 1 

ean eonstruet a Markov ehain sueh that a.s. for all n sueh that < 9^ and the 

marginal distribution of is a Markov ehain with transition probabilities 


I + 2A(7fitnessO-ir for j 2 = jl P 1, 

lP[-^n+l = 72 I^ 7l] = ^ — 2AC'fitnessO‘_ft' for j 2 = jl — 1, 

0 else. 


(7.68) 


We define a eontinuous time proeess, Z\ assoeiate to Z^. To do this, we define first 
the sequence of jump times, by Tg = 0 and 




^.7-1 “ 


— J J 


W! 


else, 


(7.69) 





























39 


where W- are exponential random variables with mean {\K{i+\){e/2)'\ {h+d+c{Ah / c))) We 

setZ^* = G [T*,T*^^). Obverse that we obtain by eonstruetion 


■ n+l 

for all t sueh that 6^^^ -(^^ 2 ) +'t < 9f. Next we want to show that 


iK 

mut. size i(e/2) 


+tJ^ 


P 


inf > 0 : Z; > \K{z + ^)(e/2)l I > = 1 - o(a,^). 


(7.70) 


Therefore, let i?f = inf{? 7 , > 0 : = \K{i + ^){e/2)~\}. Weeaneonstruet (Xj)j>i asequenee 

of independent, exponential random variables with parameter xf = \K(i + f )(e/2)] (b + d + 
c(46/c)) sueh that 

foralll<j<5f. (7.71) 

Our next goal is to find a barrier, sueh that Bf is smaller than Ui only with very small 
probability. Sinee the transition probabilities of Z® do not depend on the present state, Z® ^ — Zq 

is stoehastieally equivalent to Yl,k=i ’ where {Vi^)k&n are i.i.d. random variables taking values 
±1 with probabilities 

Wl = 1] = i + ‘^ACntn.ssCTK and P^ = - 1 ] = f - 2ACn,nes.cTK. (7.72) 
Note that E [ 14 ®] = k and |14®| = 1. Furthermore, we get 


P [Bf <ni] =F 


3f(e/4)A1 <J<nr.Y. K > f(t/4)A1 

k=l 

Hoeffding’s inequality implies that, for j > \{e/A)Kl 


(7.73) 


P 


Vk > 4ACfitnessO-Kj + i 


V2+a/2 


k=l 


< 2exp(-j“). 


(7.74) 


We take n* = eK{^ACfnnts^aK) ^ and get for all [(e/4)7f’| < j < rti, 

44lCfitnessa,,j + < r(e/4)iTl, (7.75) 

sinee 7 f“ 5 +“ ax- Then, the probability that Bf < eK{8ACfitnessO'K)~^ is bounded from 
above by 2 exp(—Tf"). Therefore, the left hand side of equation (1.10) is larger than 


P 


eir (8ACfitness'3'/c) 


Z^7 = l 


X] > Sk 


-2exp(-7f“ 


By applying the exponential Chebyehev inequality we get, similarly as in (a). 



■ err(8ACfltnessO-x) ® 


e-fi"{ 8 AOfitness'Ti^ ) 

p 

E ^1 S Sk 

= P 

1 

M 

IV 

1 


i=i 


L i=i J 


(7.76) 


(7.77) 


< exp(ii:“^;^)E 

< exp(7f“S'j^) exp fe7f(8y4CfitnessO-ir)"Mn 


xf+K° 


< exp(7f“S'j^ — eiT( 8 ACfitnessO‘K) ^CK ^’'’“)), for some small C > 0 , 

< exp(-7f“). 

This proves that P[inf{f > 0 : Z^® > \K{i + 4)(e/2)]} > Sk\ > 1 — 3exp(—iT“), and 
therefore (b) and (a) for i, provided that the lemma holds for i — 1. □ 
















40 


Proof, of (c) for i by assuming that the lemma holds for i — 1. Note that the random elements 
T*, X* and F*are not the ones of the last proof. As in (a) we couple 1) with a discrete 
time Markov Chain. Therefore, let 


Xl = \K{DtA)-\m^n))K]\ (7.78) 

and T* = and (T^)a:>i be the sequences of the jump times of {9t, 1) after 

Then, let be the associated discrete time process which records the values that XI takes after 
timeCt. sizei(./2)- 

Claim: There exists a constant C^fertative that for all j < \eK] and K large enough, 


=3 + 111'.;=i.r„+i <ef] < 2 - p+ui 

Moreover, we can choose = supjg;^ \^ - 9ic{x, x) - d 2 c{x, x)\. 

From (a) we know that the left hand side of ^ .19) is smaller or equals 


c{R^,R^) 

3b{R^) 


7-^ 86(R^) 


(7.79) 


^ - d,c{R^, R^) - d^ciR^, R^) + 0(4). (7.80) 


This proves the Claim. Note that p^(j) depends on j. Since we can choose M > (84eri™tive)(^)’ 
continuing as in Lemma [O] implies that (c) is true for i, provided that the lemma holds for 

i — 1. □ 


Proof, of (d) for i by assuming that the lemma holds for i — 1. Again we couple for 
t > with a discrete time Markov chain. Therefore, let T* = 4ar0(if) and (T^)fc>i be 

the sequences of the jump times of after 4ar(/>(ii)- Then, let (4)n>o be the discrete time 
process which records the values that OJtj ^, i.e. 


Iff Iff 

Y‘ = e 1 A’( 2 - 4) - 2 + i) + P'*') 

and 4 = 971^1(4). Define 

9f = inf > 4ar<>(i|) : 1(^0 1) - 0(*(e/2))| > Meax] A 4succ. mut. A tsity (7-82) 

Note that this 9f differs only a bit from the one defined in (b). From the proof of (b), we 
know that the density of the mutant trait has the tendency to increase. More precisely, since 
i < Qross(2/e), we have, for all - {^K] < j < , for K large enough and for e small 

enough. 


P 


= iqK] + j +1 


Y' 




K 


\ 1 I inf^sA; difix,x) 

- 2 2b 


(7.83) 


By Continuing in a similar way as in (b) with bounding the random variables in the in the other 
direction (as in (a)), implies that (d) is true for i, provided that the lemma holds for i — 1. □ 

□ 


7.3. Step 3. Similarly as in Step 2 we define a function which allows us to approximate the 
total mass of the population for a given density of the resident trait. 


Notation. Let us define 




z{R^) + (TKh{z{R^) — x) 


(r'{R^) I 

r(RX) c(RK,Rk] 


(7.84) 
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Note that 0(i/) = ^ {(f){y) -y) + 0{a\). Therefore and sinee | {vqk , 1) - (/>(Q,oss) I < 

mut. size Ocross 

Meax with probability 1 — o{aK), we get that at time size the density of the resident 
population belongs to an interval eentered at (j){C^^oss) ~ Qross '''ith diameter 2(M + [S/a] )eaK 
with probability 1 — o{aK) and henee 

, )K-^) = + 0(e4) = (/.(CLss) + 0(4) (7.85) 

mut. size C^j-oss 

with probability 1 — o(cri^). Thus, the total mass of the population also belongs to an interval 
eentered at with diameter 2{MeaK + 0(4)) < 2(M + l)€aK- 


and 


(7.86) 


Notation. Let us define 

cL. = rwc;„„)-c‘„..-e) 2 /£i(€/ 2 ) 

> Cistecs™ : Kf'i.l)- |)l < (W3)e<^/f}- (7.87) 

Note that the term — e in the definition of ensures that resident population is larger than 
^cross at time 6^^^ 

. size C^^ross ' 

First, we need a lemma for the interface between Step 2 and Step 3. 

Lemma 7.5. Fix e > 0. Suppose that the assumptions ofTheorem \7.1\ hold. Then, there exists a 
constant M > 0 (independent of e and K) such that, 

(a) Soon after size population size is close to fiC^ross ~ §)•' 


lim a j}¥ 

K^oo ^ 


^neartPiC^,,,,,-^) > ^mut. size + *S'x A 6*2 


A inf \ t > 4, = [(04. ± 3e/4)4 


K .qK 

succ. mut. ' ' ^diversity 


(7.88) 


= 0 . 




(b) A change of order efor the resident density takes more than o{a^ 


time: 


lim (TrfV 


inf * > C, ■ ®“(n,) = r± 3€/4)ii'l 


^ -L Q A a - A A Qt< 

^ ^mut. sizeCp,^,.,. ' '' ^near ip{Cg„,.,,—y) ' ^2 succ. mut. ' ' ^diversity 


(7.89) 


= 0 . 


(c) At the time when the resident density has changed of order e the total population size is still 


close to 


lim P 

K^oo 


inf <1 f > 4 : m, 1) - 


§)| >MeaK \ < d 


cross 2 

A A inf |f > sizecp,„„ ■ = [(04, ± e)K] 


(d) A change of order efor the resident density takes no more than (iax) 


-l-a/2 


■tK 

'2 succ. mut. 


= 0 . 


time: 


(7.90) 


lim (Ji^^ P 

ir-s-oo 


0^ ■ > 0mut size ^ A 0 

res. size e size / j 


K A nK ('7 911 

2 succ. mut. ' ' ^diversity 


= 0 . 


□ 


A inf : |(Ft, 1) - V'(04, - f)| > Meax 

Proof Apply the methods of of (a) to (d) from Lemma [74} 

Next, we have the following similar lemmata as in Step 2, for them let us define 

> C.sizei(./ 2 ) : - f{i{e/2))\ < {M/3)eax}. (7.92) 
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Lemma 7.6. Suppose that the assumptions of Theorem \7.1\ hold. Then, there exists a constant 
M > 0 (independent of e, K and i) such that, for all e > 0 and for all - e) (2/e) > i > 2, 

(a) Soon after 6^^ sizei{ti 2 y population size is close to 


lim arf F 

ie-s-oo 


nK ^ qK -I- a f}^ A 0^ 

^nearil){i^) ^res. sizei(€/2) ' ' ' ^2succ. mut. ' ' ^diversity 

A inf {t > ,i,ei(e/ 2 ) ■ = r(i ± l)ie/2)K] } 

(b) A change of order efor the resident density takes more than o(af^) time: 

'inf {t > C size^ie/ 2 ) ■ = [(^ ± \){e/2)K^ ] 


(7.93) 


= 0 . 


lim (tA P 

K^oo 


(7.94) 


^ ^res. sizei{el2) A dnear A (^2succ. mut. A ^diversity 


■tK 


K 


■tK 


= 0 . 


(c) At the time when the resident density has changed of order e the total population size is still 
close to A): 


lim aj} F 

K^oo 


inf <1 f 1) - ^(^(e/2))| > Meax !> < mut. (7-95) 


A A inf {t > sizei{e/2) ■ 9?l°(i>t) = [(* ± l)(e/2)/s:i } 

(d) A change of order efor the resident density takes no more than time: 

C. size (i-l)(6/2) > ...... 

A inf |f > : \{bt, 1) - f{i{e/2))\ > Mea^fj = 0. 

Proof Apply the methods of of (a) to (d) from Lemma [74} 


lim ad F 

K^oc ^ 


= 0 . 


iK A nK rq g^'j 

2 succ. mut. ^diversity 


□ 


Remark 8. Lemma [73] and [7)^ imply that the density of the resident trait deereases to the value 
e. Moreover, 


P 


iK 

res. size e 


> size A A 0lersityl = (7.97) 


and 


P 


UgK , 1) - f(e)\ > Meax = o(aK)- 


res. size e 


(7.98) 


7.4. Step 4. After the time 0l we have to wait less than ln(A')crl“^^ time to know that 
the resident trait is extinet with high probability. 

Notation. Define = inf{f > 9g : \ {i>t, 1) - ^/(0)| < (M/3)eai^}. 


Lemma 7.7. Suppose that the assumptions of Theorem \7.1\ hold. Then, there exists a constant 
M > 0 (independent of e and K) such that, for all e > 0 

(a) Soon after 0l the total population size is close to '^(0).' 


lim ad F 

K^oa ^ 


nK ^ nK -U A 0^ A 0^ 

^nearijj{0) ^res. size e ' ^K ^2 succ. mut. ^diversity 

A inf {t > 01 gg) : = [(1 ± |)eiPl } 


(7.99) 


= 0 . 
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(b) A change of order efor the resident density takes more than o{aj^ ) time: 


lim aJ 

K^oo ^ 


P 


inf {t > ef,, = [(1 ± l)eK ]} 


^ ^res. size e “ 1 “ A Onear'{p{ 0 ) A 62 


K .oK 

succ. mut. ' ' ^diversity 


= 0 . 


(7.100) 


□ 


Proof. See proof of Lemma |7.4| 

Lemma 7.8. Suppose that the assumptions of Theorem \7.1\ hold. Then, there exists a constant 
M > 0 (independent of e and K) such that, for all e > 0 


lim (Tj} V 

K^oo 


qK 

^ res. size 0 


K .nK 

2 succ. mut. ' ' ^diversity 


> (»Lrm + HK)<yK^°A A » 

A inf {t > ■■ \{vu 1) - i/'(0))| > Milok} 


(7.101) 


= 0 . 


Proof. To prove this l emm a we use a eoupling with an eontinuous time branehing proeess as 
in the proof of lemma For any < t < ^f^ucc. mut. A tsity A inf{t > : 

I (Ft, 1) — -0(0)) I > Meax}, any individual of lHT°(Ft) gives birth to a new individual with trait 
with rate 


(1 - UK m{R^))b{R^) e [b{R^) -uxb, b{R ^)], (7.102) 

and dies with rate 

d{R^) + c{R^, R^)m%i)t) + [ c{R^,0dM0, (7.103) 

JxxN 

whieh is larger than dz = d{R^) + c{R^, R^ + axh)z{R^ + cxh) — where 

^Maideath = M + c[3/q;] — 2 hd2c{R^, R^). Therefore, we eonstruet, by using a standard 
eoupling argument, a proeess Zt sueh that 

> (H°(Ft) (7.104) 

forall0,^^^^(o) < t < Cucc.mut. ACersityAinf{f > ^ I 1) “ ^(0)) | > Mcax}. The 

proeess Zt is a linear birth and death proeess starting at [|eiF], with birth rate per individual 
bz = b{R^) and with death rate per individual dz- Sinee 

bz-dz = fiR^,R^ + axh) + C'.^.ideath^cr^ (7-105) 

= -axhdif{R^ + axh, R^ + axh) + C'j^aideath+ 0{{axh)‘^) = -ax^K 


is negative and of order ax, the proeess Zt is sub-eritieal. Note that ^x > > 0. 

Let rf be the first hitting time of level i hy Zt, then we have 


^[^f 2 eK] < ^o] < exp(-iT“) (7.106) 

eompare with the proof of Proposition |9.5[ Sinee Zt > QJl°(Fi), we obtain also that, with high 
probability, 91T°(Fi) stays smaller than \2eK~\ before it dies out. For any t > 0 and n G N, the 
distribution of the extinetion time of Zt for bz f dz is given by: 


/ dz - dz exp{{dz - bz)t) \ 
\bz - dzexp(((iz - bz)f) ) 


IPn(r(f < t) 


(7.107) 
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(cf. [l2l P- 109 and [0). Therefore, we ean eompute in our ease where dz — bz = (Tk^k with 
uniformly positive 


P 


To <ln(it:)cT^ 


-l-a/2 


dz - dz exp {{dz - bz) \n{K)ap^ j 
bz - dzexp {{dz - bz) \n{K)a~[^~°‘^ 


, , /- —cx.j 2 

dz — dzK^^^ ^ 
dz — ck^k — dzK^^'^K 


ieK 


= 1 - 




^eK 


dz{K^^'^K — 1 ) + cr k^k , 


> {l-aK{\eK)-^K-Y^^^ 

> 1-0{(TkK~^)>1-o{(Tk), 




(7.108) 


whieh proves the lemma. 


□ 


7.5. steps. After the extinction time of the resident trait, we have to wait at most ln(7f)cr^^ 
time until the population is monomorphic with trait + axh. 

Lemma 7.9. Suppose that the assumptions of Theorem \7.1\ hold. Then, there exists a constant 
M > 0 (independent of e and K) such that, for all e > 0 


lim (Jj} ¥ 

K^oo ^ 


iK 


^fixation ^ i^res. size 0 ^ ^2 race. inut. 


A 

' ' ^diversity 


-l-a/2- 
K 

A inf {t > : |(i>i,ll) - f{0)\ > Meax] 


(7.109) 


= 0 . 


qK 


|Supp(i/f)| = - 


Proof By the last lemmata, we have 6'^ation = inf{f > 
i/’(0)| < {M/‘f)eax} with probability 1 — o{ax)- Set D = {k E N : 1 < iM’^{i>QK ) < 
eaxK}. Then \D\ < [3/a] and none of these traits are successful since we have seen that 
^res. size 0 Smaller than and with a probability of order 1 — o(cTii'). By applying 

Proposition |9.3| and using the Markov inequality, we obtain that the life time of each of these 
subpopulations is with probability 1 — o{ax) smaller than \n{K)af^~°'^‘^. Therefore, if no new 


mutant is bom betw een 6 ^,^ g and 6^^ Q+\n{K)a^ we obtain the claim. On the other 


hand, as in Lemma 6.4 the number of mutants born in the time interval [6^^,. 9^^ ^j^gg + 

ln(7f)cr^^~“^^] is stochastically dominated by a Poisson point process, A^{t), with parameter 
a uxK, where a = sup3,g_:^. z{x)b{x)m{x) + 1. Hence, the probability to have no new mutant 
in this interval is 


P[A^(ln(A:) 


cr 


-l-o/2^ 
K 


) = O] = exp(— \n.{K)a^ uxK) > exp(—a^^^) > 1 — o(l) 


q/2n 


(7.110) 
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Because the probability that a mutant is successful is of order ax, the probability that a suc¬ 
cessful mutant is bom between times 6^^ ^res. sizeo + o^ax)- Since 

(7.111) 


F[A^{\n{K) 


a 


-l-a/2. 


K 


< rs/al] 


exp (— ln(7f)cr^^ auxK^ 


[3/01 , , -l-a/2 

m(A)cT^ auKJ^ 


j=0 


> 1 - 


> 1 - = 1 - 0{(Tk), 


there are maximal [3/a] unsuccessful mutations in this interval. With the same argument as 
before the life time of each of these subpopulations is with probability 1 — o(ai^) smaller than 
^-i-a/4 „,;^v 1 _ o{a k) the maximal possible time interval where 

-l-ajA 


\n{K)(T^‘^ Therefore, with probability 1 


at least one mutant individual is alive is s mall er or equal \n{K)aj^ -f [3/a] hi{K)a~p^ 


<C \i\{K)a 


- 1 - q ;/2 
K 


Recall from Lemma 


7.7 


that if 1) — V’(0)l > {M/3)eaK at the first 


time when the population is again monomorphic, then the time the process needs to enter 
the (M/3)e(Tii'-neighborhood of '0(0) is smaller than Sk, which can be chosen smaller than 
(T^“/ (Kuk)- This proves the lemma. □ 


This ends up Step 5 and the second invasion phase. Note that the estimates of the two phases 
do not depend on the exact trait value of the resident trait, especially the a priori different 
constants M. In fact, we can use in all lemmata the same constant M, namely the largest. 
Therefore, we can apply our results for the successful mutant trait Rf = + axh, which 

is the next resident trait by using the strong Markov property for (z/, L) at the stopping time 

nK 

^fixation* 


8. Convergence to the CEAD 


Our goal is to find Tq > 0 and to construct, for all e > 0, two measure valued processes, 
(/r]’'^’^, f > 0) and (/rt’'^’0 1 > 0), in D([0, cxd), A'l(A’)) such that 


lim P 

K^oo 


Vf < 


Tn 


— Ku 


k<^k 


and forj e {1, 2} 


l,K,e j K j 2,N,e 


= 1 , 


( 8 . 1 ) 


lim P 

iT—>oo 


sup 

0<t<To 




>6{e) 


= 0 ., 


( 8 . 2 ) 


for some function 5 independent of x,K such that ()(e) —)■ 0 when e —)■ 0. This easily im¬ 


plies ( |4.5[ ) for all T <Tq. 

The result for all T > 0 then follows from the strong Markov property. Indeed, the construc¬ 
tion below implies that there exists a stopping time r G [Tq/2KukO'‘^k, Tq/Kukc^W (a fixation 
time) such that, with probability converging to has a unique (random) point Y as support 
and a total mass belonging to [z{Y) — Max, z{Y) + Max]- Hence ( |8.1[ ) and ( |8.2[ ) also hold for 
the process (^'/[,_t, t > 0), and ( |4.5[ ) is thus true for all T < 3To/2. We obtain (4AI for any fixed 
T > 0 by induction. 
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8.1. Construction of two processes and such that ^ ^ Fix T > 0. 

Let denote the random time of i-th invasion (i.e. 9f = 6*j^nvasion)’ ^i^xation the time of i-th 
fixation and Rf the trait of the Ath sueeessful mutant. Let us fix the following initial eonditions 




_ jdK 

Kq — Kq 


Aax, Rq’^ = Rq +A(Tk and = 9^’^ = 0. Assume that we have eonstructed 


nir,i 


nir,2 


and 9f^’^, and R. 


K,1 


K. 2 

and ’ . By Theorem 6.2 and Markov’s property, we ean eonstruet 


two random variables Rf^ and Rf^ sueh that 




K,1 


■i+l 


p^il ^ tdK 


jdK ^ p^>2 


R. 


.K,2 


with probability 1—o(crii'). Moreover, Rfj^\- 
1 - 0(e) and Rf’l - Rf'\ ^ 

Corollary |6.10| ) 


p^a_ pir 

-iij — 

< A(Tk- The distributions of R, 


R^ = i?: 


K,2 


.K,l 

i+l 


i+l 
K+ 


(8.3) 

-Rf'^ with probability 


Rf'^ and Rff. - are (ef. 


K,2 


r\{Rf,h)^¥[R^;\ = Rf+akh] = 


M{Rf ,l)ql{Rf ,1) ^ ^ _ p\{Rf) 


and 


rl{Rf,h)^nRt;i = Rf+<ykh] = 


- rorD^.2 


K 


pliRf) 

M{Rf,h)qliRf,h) 

M{R^,h)ql{R^,h) 

PliRf) 

M{R^ ,A)ql{Rf ,A) 

PURf) 


Pl(Rf) 


\+l 


if h = 1 


if he {2,..., A} 


(8.4) 


+ 1 - 


pliRf) 

pliRf) 


ifh e {1,..., A—1} 
ifh = A 


where 


q\{x, h) = h 


dif{x,x) 

b{x) 


a 


Bernoulli^) 


b[x) 


2 , 
Bernoulli^ 


(8.5) 


( 8 . 6 ) 


and p^j{x) = Qji^^ h)M{x, h) for j = 1, 2. (Note that we ehanged a bit the notations 

of Corollary 6.10 to make explieit the dependenee on e and Rf.) Sinee we assumed that the 
fitness gradient dif{x, x) is positive and uniformly lower bounded on X, the transition proba¬ 
bilities r^(a;, h), j = 1,2 are uniformly Lipsehitz-eontinuous funetions of x with some Lipsehitz 
eonstant OLip. 

By Theorem 6.2 and Lemmata 6.7 and 6.4[ we ean eonstruet two exponential random vari¬ 
ables, eI^I and E-_^l ,with parameters af’'^{Rf)pl{Rf)aKUKK and a 2 ''^{Rf)p\{Rf)<JKUKK 
given by 


(h''^{x) = {z{x) — ea KM)b{x)m{x) 


K,t 


x) = {z{x) -f eaxiM + |'3/a]))(6(a:)m(a:) -f Aax 


sueh that 


P 


p^s',2 ^ nK 

^i+l — 9 i+l 


affixation < E^+l + ln(iT)a-^-“/^) = 1 - 0{aK). 


(8.7) 

( 8 . 8 ) 

(8.9) 


Note that this inequality involves 6*j^xation instead of 9f- sinee we apply the Markov property at 
the fixation time of Lemma |7^ before we ean apply Theorem |6.2[ However, Lemma [7^ entails 
that we also have 


P 


Ef+l < Of+i 


< Eli; + 6ln(K+X‘‘'") = 1 - o+kI 


We then define 


gKl _ gK., _ 


= bS + 61ii(A')<t- 


-l-a/2 


and 


oRfi _/3^’2 _ jpRP 

9i+i — -t^i+i- 


( 8 . 10 ) 


( 8 . 11 ) 
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In addition, by their construction in Section m it is clear that the random vectors 

— p^’^w ofo { pK\ 


\+l 


\+l 


Rf' )}i>o are independent conditionally on {Rf) 


■f )i>o- 


Lemma 8.1. With the previous notations, the stochastic processes, and , in 

D([0, oo), M.{X)) defined for allt > 0 by 


p]’^ = {z{Rf) - {Me + C)a^) Vi, 


/ortG[0f’\Ci)n[0f,0f+,), (8.12) 
Pt’^ = {z{Rf) + {Me + [3/a]e + C)aK)5^Ka, fort e [^f’^ O^fl) n [Of, (8.13) 

for some constant C independent of K, x, e, satisfy for allT > 0 


lim P 

K^oo 


Vi < 




1,K j K , 2,K 

pf ^ dt 


= 1 . 


(8.14) 


Note that the support of p^’^, j = 1, 2, is defined from the sequences {Rp’^)i>Q and {Op’^) 
but the mass of p^’^ is defined from the sequences {Rf)i>o and {9f)i>i. 


i>l 


Proof Let us fix T > 0 and F > 0. Since each of the steps previously described holds with 
probability 1 — o{aK), we deduce that the above construction can be done on a so-called good 
event of probability 1 — o(l) for all integers i < T/ap- Since in addition {x)p 2 {x) is uni- 
formly lower bounded by a positive constant a on A’, the random variables E- ’ can be coupled 
with i.i.d. exponential ones of parameter oKukctk, and hence P[ < T/{Kupcrf) ] 

is smaller than the probability that a Poisson process with parameter aKupCK is larger that 
[F/cr^J at time T / {Kupcrj^)- By the law of large numbers for Poisson processes, we deduce 
that, provided that F > Ta (which we assume true in the sequel). 


lim P 

K^oo 


nK,2 


< 


T 

Kupaj^ 


0 . 


(8.15) 


Let us recall that, on the previous good event of probability 1 — o(l), the number, the trait and 
the size of the living mutant populations and the size of the resident population are controlled at 
any time in the Lth first phase (Lemmata |6.3| and 6.91. In addition, during the Lth second phase, 
the number, trait and size of living mutant populations are controlled (see all the Lemmas of 
Section 1^, the total mass of the populatio n sta ys withi n the Meap-neighborhood of (j){y) or 


fi{y) for some y G [0,z{Rf)] (Lemmata 7.4 and 7.6). Since |0(r/) — z{Rf)\ < Cap and 
\'4>{y) — z{Rf)\ < Cap for some constant C, as seen in ( |7.31[ ) and ( |7.84| ), and since the 
sequences {Rf^)i>o for j = 1,2 and {Rf)i>o are all increasing on the good event, we deduce 


the required comparison between the supports of p^ 
event. Since we used z{R^, 


r,K p 
of 


and p^p^ for t < 


if) to define the masses of p]'^ 
between the masses is also clear. 


and Pf- 


2,K 


, the required comparison 

□ 


Note that, since the function ^ may not be non-decreasing, replacing z{Rf) by z{Rf'f in 
the definition of p]’^ may not imply the required comparison between the masses of p\'^, of 
and 

The next goal is now to prove the convergence of both processes pff 2 for j = 1,2 to 
z{xt)Sxt in probability in L°°{M{X), || • ||o). For this, we will use standard convergence results 
of Markov jump processes. However, the two processes p^’f j = 1,2 are not Markov because 
the Lth jump rates and transition probabilities defined above depend on Rf which is close, but 
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different from Rf'^. Therefore, we introduee a small parameter rj > 0, and we shall eonstruet 
two Markov proeesses j = 1,2 in D([0, oo), A^(A’)) sueh that 

Z. _A 

KuKal " 


lim P 

K^+oo 


j 1,K j K , 2,iC . K^2,€^rj ^ 

^^{t-l/{KUKCTK)RO ^ ^ ^ 5 Vf < 


= 1 , 


(8.16) 

where is the first time where the distanee between the support of and is larger 

than T]. The last equation will be proved below in Seetion [8^ The time-shift of —1/{KukO'k) 
in is due to the terms Q\n.{K)a~Z~°‘^‘^ in ( |8.11[ ). We will next study the eonvergenee 

of these two Markov proeesses when K ^ oo and prove in Seetion 8.3 that, for a eonvenient 
ehoiee of r], there exists some Tq > 0 independent of K, x, e, rj sueh that 

To 


lim P 

K^-\-0O 


oK , 

. " KuK<yl, 


= 0 . 


(8.17) 


8.2. Proof of ( |8.16| ). For all x G A", we define (rj’’'(x, h),l < h < A) and {r 2 ^{x, h), 1 < 
h < A) by, for all 1 < £ < A, 


^ri’^(x,/i) = 


h=i 


'Y^{r\{x,h) + Cl^T]) 


h=l 


and 




Note that 


h=l 


Yl(rl(x,h) - Cl^ri) 


h=l 


Al> sup '‘^r\{y,h) (8.18) 

y&[x,x+rj\ 


V 0 < inf ^r^(|/,/;,). (8.19) 

y&[x,x+rf\ f—^ 
h=l 


UC, 


and r 




x, 


are probability distributions on {1,..., A} for all x G A and 


that, by standard eoupling arguments, for all x < ?/ such that y — x < rj, the distribution 
fl''^{x, ■) is stochastically dominated by the distribution rl{y, ■) and the distribution r2(x, •) is 
stochastically dominated by the distribution f 2 ^{y, •). We define similarly 

K,e, 


— K^e^ri/ \ K^e 


1 (a;) = fli - Cli V < inf oi ’\y)pl{y), 

^ y&[x,x+y\r\X 


( 8 . 20 ) 


and 


df’"’''(x) = af’"(x)p 2 (a:) + (PLip// > sup a^’"(|/)p 2 ( 2 /), 

y&[x—ri,x\r\X 

where C^ip is a uniform Lipschitz constant for the functions j = 1,2. Note that 


K,ei 


( 8 . 21 ) 


a 


K,e,7]. 


1 


x) > 0 for all X G A if r; is small enough. 


It is then clear that there exist two Markov chains (R, 


K,j,v\ 


}i>0 


, j = 1,2, with initial condition 


T3K,j,y _ 

/to 


- 

— /to 


and with transition probabilities rf’'^’^{x, h) from x io x + h, such that, for all 


i > 0 satisfying i?- 


K,2,r] 


Sir,lx ^ 


S i?, 




K,l,r] 




_ pTk _ 

2+1 aiivi 

Similarly, there exists random variables E, 

conditionally on {R, 


■R. 


K,1 


and R, 


K,2 


pir,2 ^ TjK,2,ri 
+ ^ ++1 


jyK,2,ri ^ 

■ Ttj ^ ++1 


Rf'\ (8.22) 


tributed with parameters 


ZZ, j = 1,2, independent and exponentially dis- 


such that - 


< E, 


K,2 

i+l 


and E^l < Ej^^. We then define 


9. 


K,j,ri 


= E., 


K,j,e 

i+l 


with 6*, 


K,j,ri 


= 0 . 


Since the function z is (^Lip-Lipschitz, it is clear that (|8.16) is satisfied for the processes 


■^Lip 

K,l,e,y _ vK,l,y 


and 


Pt 

X,2,e,r7 

Pt 


= iz{Xr’^) - (Me + C)aK - 


(8.23) 


— ’ ’^) + (Me + [3/Q;]e + C)(Jk A ^ (8.24) 
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where 

^K,i,v ^ foj. t ^ ^QK,i,r, ^ + Q{^ + 1) ln(ir)a“^-“/'), (8.25) 


and 


X, 


K,2,t) _ ^K,2,r] 


for^G[0T’^’^0S’")• 


(8.26) 


By construction, the processes and are Markov jump processes, but the process 

XK,i,r] jg because of the terms 6 ln(X)(j^^~“^^ involved in its definition. However, the 
process = (^(Xf’’^’’') — eaxM — Cl^pr])5^K,i,r, is Markov, where 


j^K,1,7] _ J^K,l,ri 


fortG[0T’^^0S’^). 


(8.27) 


The proof of (|8.15[) above also applies to the processes and Since in addition 


the support of is non-decreasing, it follows that /r 

t <T/{KuKcrj^) with probability 1 + o(l). Our assumption ( |4.2| ) entails ( |8.16[ ) 


-2-a/2 


K,l,e,7] 

(i-er \n{K)a 


)V0 




A',l,e,?7 


for all 


8.3. Convergence of when K —)■ +oo and proof of ( |8.17[ ). The two Markov processes 
^i^/’{Kuk^k) ^i^/’(Kuk<tk) to the framework and assumptions of Theorem 2.1 of 

Chapter 11 of Ifl^ : their state spaces are (up to a translation) a subset of and their transi¬ 
tion rates from 2; to z + hax have the form a]^[l3h{z) + 0{aK)] for some Lipschitz functions /S/j. 
For such a process X, provided Xq converges a.s. to xq, the process > 0) converges 

when cT/i' —)• 0 almost surely in L°°([0, T]) for all T > 0 to the unique deterministic solution of 
the ODE dx{t)/dt = h/3h{x) with a:(0) = xq. In our situation, we obtain, for j = 1, 2, that 


lim sup 

K ^+00 tg[o,T] 


X 


K,j,rj 

t/iKuKcrj^) 


Xj{t) 


0 a.s.. 


(8.28) 


where Xi and X 2 are the unique solutions such that Xi(0) = ^2(0) = x of the ODEs 

= [z{xi{t))b{xi{t))m{xi{t))pl{xi{t)) - ^ /ifi’^(a:i(f), h) (8.29) 


and 


dx 2 {t) 

dt 


[z{x 2 {t))b{x 2 {t))m{x 2 {t))pl{x 2 {t)) + Cupp] h). 


h=l 


Lemma 8.2. For all T > 0, and for j = 1,2, 


sup \xj{f) — Xt\ < CTe^ (?7 -f e), 

i6[0,T] 


(8.30) 


(8.31) 


for a constant C independent of x, T, e and p, where Xt is the solution of the CEAD ( |4.3[ ) with 
initial condition xq = x. 


Proof We only write the proof for j = 1, the case j = 2 being similar. Since the functions 
j = 1,2, b, m and pi are bounded by constants independent of X, e, p, we have for all 
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t E [0, T] and for a constant C > 0 that may change from line to line, 


\xt-xi{t)\ < + 


{zbmpl){xi{s)) ^ hrl’'^{xi{s), h) 


h=l 


- {zbmp\){xs) ^ 


h‘^M{xs, h)dif{xs,Xs 


h=l 


b{xs)pl{xs 


ds 


< + ACl^)Tr^ + C \xs- x,is)\ds 


nt ^ 

h=i 


rl{xs,h) 


hM{xs, h)dif{xs,Xs) 


b{xs)pl{xs) 


ds, 


(8.32) 


where the last inequality follows from the uniform Lipsehitz-eontinuity of all funetions involved 
in the eomputation. Now, \pl{x) — pl{x)\ < Ce and Pj{x) > c > 0 for j = 1,2, for some 
eonstants C,c> 0 independent of e and x. Henee, there exists a constant C such that 


\xt — Xi{t)\ < CT{p + e) + C / |x5 — Xi(s)|(is 


/•i A 
>^0 h=l 


q{{xs, h)-h 


dif{Xs,Xs 


b{xs 


M{xs, h)ds. (8.33) 


In view of ( |8.6[ ), we obtain \xt — a^i(f)| < CT{rj + e) + C f* \xs — Xi{s)\ds. Gronwall’s lemma 
ends the proof of Lemma [8^ □ 

In view of Lemma [8^ there exists Tq > 0 independent of x,e,p sueh that, for all ?7 > e, 
supigfo.Toi ~ ^t\ < v/^- Let us fix ?7 = e. Combining ( |8.28| ) with the last inequality 
entails ( 8'l7| ). 

8.4. End of the proof. 

Proof of Theorem l?^/] Defining ^nd and combining ( |8.16[ ) 

and we see that we have defined a eonstant Tq > 0 sueh that 


lim P 

K^-\-oq 


A''(t-i/(ir«if(T^))vo ^ Pt ^ ^ d't ^ Pt , vr 5; Kukctk 

,e 

1/(Kuk(tk))'^0 


1,K . K , 2,K , K,2,e 


Tn 


= 1, (8.34) 


This is ( 18.1[ ) with = P^t-i/(KuK(TK))vo- L only remains to eheck ( |8.2| ). 

Using that x] = e, we get 


(8.35) 


< U 


+ 


e + <JK 


< c 


z{xt) - z I + kt - - (^k) V 

Tt-Z)V0lKu^aj, - ^l((^ - ^k) V 0)|] 

+ aK+ sup (|x(t_<,^)vo - a^il + |xi - Xi(f)| + 
xe[o,T] ^ 
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for some finite eonstants C,C' > 0. The analogous estimate holds for for 2 • Setting 

for example 5(e) = ^/e, ( |8.2| ) follows from ( |8.28[ ), Lemma [8^ and the uniform eontinuity of Xt- 
This ends the proof of Theorem |4.1[ □ 


9. Appendix 


In this seetion, we state and prove several elementary results, whieh we used in the proof of 
our main theorem. Reeall that || . ||g is the Kantorovieh-Rubinstein norm on the veetor spaee of 
finite, signed measures on X, i.e. 




A’ 


fdfit : / G Lipi(A’) with sup |/(x)| 



(9.1) 


where Lipj^(A’) is the spaee of Lipsehitz eontinuous funetions from X to M. Let Adir(T:’) be the 
set of non-negative finite Borel-measures on X. 


Proposition 9.1. Let{i>^,K > 0} and fj, be random elements in H If, for all 

5 > 0 , 


lim P sup ||z/f -/iJo > ^ 

ir-s-oo |_ 0<t<T 


0 , 


(9.2) 


then converges in probability, as K ^ 00 , with respect to the Skorokhod topology on 
D([0,T],M(A)) top. 


Proof Let us equip with the topology of weak eonvergenee. Obverse that this topology 

is metrizable with the Kantorovieh-Rubinstein norm, see flU Vol. II, p. 193. Let A be the elass 
of strictly increasing, continuous mapping of [0,T] onto itself. If A G A, then A(0) = 0 and 
A(T) = T. The Skorokhod topology on D ([0,T], (AIf( 5b), || . ||o)) is generated by the distance 


d(/i, u) = inf 
AeA 


max 


sup |A(f) 

ie[o,T] 


sup \\pt 
ie[o,T] 



(9.3) 


on D ([0,T], (AfF(5b), II . ||g)), see e.g. Il3l. Chap. 3. Since the identity lies in A it is clear that 
d{p, 0 ) < supigjo.T] Wdt — Therefore, if a sequence of random elements with state space 
D([0, T], AfF(5b)) equipped with the metric induced by the norm supjgjg.T] llf'-t llo convergences 
in probability to p, it also convergences in probability to p if D([0, T], AIf(A’)) is equipped 
with the metric d. □ 


Proposition 9.2. Fix e > 0 and let ax a sequence in K with K 92 +« ctf 1- Let be a 
Markov chain with state space Nq and with the following transition probabilities 

{ 1, for i = 0 and j = 1, 

I - CiiK~^ + C 2 eaK, fori>l and j = i + 1, (9.4) 

I -f CiiK~^ — C 2 eaK, for i > 1 and j = i — 1, 

for some constants Ci > 0 and C 2 > 0. Let Ti be the first hitting time of level i by Z and let P^ 
denote the law of Z conditioned on Zq = a. Then, for all M > ^ and for all a < ^MeaxK 

Pa [riMeai^K] < T-q] = 0. 


(9.5) 
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Remark 9. The proposition can be seen as a moderate deviation result for this particular Markov 
chain. More precisely, we can prove that there exist two constants M > 0 and > 9 which 
depend only on Ci and C 2 such that for a < ^MeaxK 

Pa [r\MeaKK] < To] < cxp {{\MeaKKf - a^)) , (9.6) 

for all K large enough. 

Proof. We calculate this probability with some standard potential theory arguments. Let h\MecrKK^,fl{a) 
be the solution of the Dirichlet problem with A = 0, i.e. 

^^h^MeaKK^fi{.x) = 0 for 0 < o: < \MeaKK~\ 

h\MeaKK^fl{x) = 1 forX > 




'.x) = 9 


for X = 0. 


Therefore, we obtain for 0 < a < \MeaKK~\ (cf. [|5l p. 188) 

E d _i _ 1 

i=l 7r(2) 1) 




E l IVli 
2 = 1 


\Mec7KK^ 1 


(9.7) 


(9.8) 


7r(2) p(i,i—l) 

where tt = (7r(0), 7r(l), 7r(2),...) is an invariant measure of the one-dimensional Markov chain 
Zn- In our case any invariant measure vr has to satisfy, for alH > 1, 

7r(0) = p(l, 0)7r(l) and 7T{i) = p{i — l,i)7r{i — 1) + p{i + l,i)7T{i + 1). (9.9) 

Therefore, vr with 7r(0) = 1, 7r(l) = and 7r{i) = pO-y-lj the unique invari¬ 

ant measure for the Markov chain Zn. Thus we get from ( |9.8 ) 

-r-ri-l p{j,j-l) 

Z^i=l ll7 = l 


,0 (®) 


E l IVl ( 
2=1 


ti=i p(i,i+i) 

\MeaKK] tti-I p(j,j-l) 


ni=l p{j,j+l) 

i._ f l-\-2CiK~^j—2C2^crx 
y l—2CiK~^j-\-2C2^o-K 


E“=iexp E;=;in U^ 


E l IVl ( 

i=l 


\MtaKK] 


exp 


S^i-l 1 / I+2C1A' ^j-2C2eaK \ 

^i=l Vl-2Ciir-lj+2C2e(7K J 

"-V-' 

= -f{j) 


(9.10) 


For all j < MeaxK we can approximate / (j) as follwos 

4CiX-ij-4C2«:a-K A _ iCiK-^i-AC2e<TK 


f{j) + l-2CiK-^j+2C2eaK J l-2CiK-^j+2C2e(7ii 

= 46^1;^ — 4C'2ecri^ + O ((;^)^ + 

= AC,^-AC2eaK + 0{{MeaKY) 


O 


ACiK ^j-4C2£g-K 'i2 
l—2CiK~^j-\-2C2^o'K 


Therefore, 


E :=,^ MET = i ^ c , j ^+ omecTK ? 




ES"’"'' exp(j:‘-i 4Cii - 4C2(aK - O ((M(Uk)- 

aexp(2C'ia^7f“^ -I- O {a^MeaxY) 

I Mi 
2^1=1 


< 


(9.11) 


(9.12) 


\M<LaKK~\ 


exp {2CiK-^{i^ 


< 


i) - AC2eaKi - O {{i - l){MeaKy)) 
aexp{2Cia'^K~^ + O {a^MeaxY) 


“P i'iCiK-V - (2C,A--‘ + - o HhU^KY))' 


E l Ai e 

i=i 
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Choosing M > ^, if a < then 


aexp(2C'ia^i^ ^ + O (a(Me(Ji^)^) 

\Mea^K]fi{a) < 1 ^MeaxK] exp {{^CiM - 2C'2)MeV|.iC - O {{eaMfK + ea^ 

< 2a([MecTi^i^'])-^ exp ( 20 ^ - \{\MeaKK'\f)) 

< exp {-CsR-^ {{l\MeaKK]y - a^)) . 


Since K ok when K tends to infinity, (9.51 follows. 


M)) 

(9.13) 


Proposition 9.3. Let (Zt)t>Q be a branching process with birth rate per individual b and death 
rate per individual d. Let Ti be the first hitting time of level i by Z and let Pj denote the law of 
Z conditioned on Zq = j, and Ej the corresponding expectation. Then 


Pj [Tk < T-q] = 


Pl[T-fc < To] 


[b - d] + 
b 


< 


El [Tk A To] < 


{d/by - 1 
{d/bfi -1 

k-^ 

1 + In(fc) 
b ’ 


for all 1 < j < k — 1, (9.14) 

and (9.15) 

(9.16) 


where \b — (i]+ = max{6—d , 0}. Moreover, if Zt is slightly super-critical i.e. b = d + e, then 


En[rfc A To] 1 + In(fc) 
max —-p-r < - 

n<k P„ Tfc < To e 


(9.17) 


Proof Let pj = Pj[rfc < tq]. Then po = 0, = 1 and pj = ^ p^+i + ^ p^.i for all 

I < J < k — 1 by the Markov property. From this recursion, we obtain the characteristic 
polynomial 

P{x) = bx‘^ — {bd)xd. (9.18) 


With its roots 1 and d/b, we obtain the following general solution for the recursion 


Pn 


Kq ■ + Ki 



(9.19) 


where kq and ki are constants. From the initial condition po = 0 and pk = 1, we obtain 
kq = —{{fy — 1)”^ and Ki = {{fy — 1)“^- Therefore, 


d\n 


Pn = 


() 


ay-I 


and Pi = 


f-1 


yy-1 i + | + ... + (^)fc-i 


(9.20) 


Ifd > b, this computation implies thatpi = Pi[Tfc < tq] < 1/k and [b — d].^. = 0. If d < b. 


< To] - 


b-d f-1 , . u 


ity-i a-my 




yfy-i " ^\iy-i ay-I 
f(i-l) 1 1 


1 + --- + 0 


b\k 



(9.21) 
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Similarly, if Cn = E„[rfc Aro], then is the solution of the following non-homogenous Dirichlet 
problem: 


en = —1, for n G {1, /c — 1} 

Cn = 0, forn e No \ {1, - 1}, (9.22) 

where {^f){x) = x(b[f{x + l) — f{x)]+d[f{x — l) — f{x)]) is the generator of the branching 
process Z. Therefore, we have to solve the following non-homogeneous recurrence 

en +2 - ^e„+i + fen = and cq = = 0 (9.23) 

We solve this by variation of parameters. Thus, we first solve the associated linear homogeneous 
recurrence relation: 

hn+2 - + ihn = 0 (9.24) 

As we have seen before /r„ = ^ 2 1 + K 3 {fy for any ^ 2 , ks G M solves the equation. Obverse that 
this functions are the harmonic functions of Second, we have to find a particular solution. 
Let {xij, X 2 j) the solution of the system of linear equations 

x„ + (iy*'xy = {) (9.25) 

= -kAj. 

then 


e 


p 

n 


n—1 


i=o 


n—1 

J=0 


d\'n' 
b. 




(9.27) 


is a particular solution. Now, we obtain we obtain the following general solution for the recur¬ 
rence: 

Cn = /in + = ^2 + K3(f)” + ^ ^ - ij. (9.28) 

i=i 

We have the boundary condition eo = = 0, therefore K 2 and ^3 are given by the solution of 

the following system of linear equations 



(9.29) 

(9.30) 


and we obtain that 


1 





- 1 


1 



(9.31) 
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With this formula we can easily prove the second inequality of the proposition, 


ei = 




n=l ^ 0' 


n=l 


1 + ln(/c) 


(9.32) 


Finally, we obtain for slightly super-critical Zt, i.e. with b = d + e. 


^n[rk A To] _ ^ 
^n[rk < ro] Pr. 


-Y- 

b—d / j A 


i=i 


d 


k-j 




1 iHY-m-im 


<1 






<0 


< 1 1 < 1 

~ e ^ j ~ e 
i=i •’ 


which proves (9.11). 


(9.33) 

□ 


Proposition 9.4. Let {Z^)t>Q be a sequence branching process with birth rate per individual 
6 > 0 and death rate per individual d > 0 and \b — d\ = 0{aK), where ax 

Let Ti be the first hitting time of level i by Z and let denote the law of Z conditioned on 
Zo = j. 

(a) The invasion probability can be approximated up to a error of order exp(—iF"), i.e. 

lim exp(if")|Pi < ro] - LE 

K^oo 0 


= 0 . 


(b) Ifh > d (super-critical case), we have exponential tails, i.e. 

-Q:/3\Tn, r -l-o/2 


lim exp(a^"'^ )Pi YctkK^ >\'n(K)ap^ 

K—¥0Q 


and 


Proof, (a) Compare with ( |9.14[ ) that 

Pi [r\eaKK^ < T-q] = 


^ 'fo 

lim exp{K^)F^,,^K] Yk] > To] = 0 

(d/b) - 1 


= 0 


K^oo 


(9.34) 

(9.35) 

(9.36) 

(9.37) 


{d/b)^'^^KK] _ i' 

If b > d (sub-critical case), there exist two constants C™'’ > 0 an d 67™'’ > 0 such that 1 -f 
< d/b <1-1- 67™’’ (Tx- Therefore, the left hand site of (9.37) does not exceed 


< 


67™^^ 


= o(exp(-iT“)). (9.38) 


(1 + C™'’aK)r-Ki^l _ 1 - exp(C™'’ai^ [eaxiTl - 67(a|eiT)) - 1 

The last equality holds, since K‘^°‘ <C cr^-iT. If 6 > d (super-critical case), we obtain similarly 


Pl[T-fc < To] - 


(b) Compare with [HI page 41, that 


b-d 


f-1 

b 




= o(exp(-iT")). 


(9.39) 


Pi 




'f\eaKK~\ ^ Fl 


(9.40) 


< exp — 


ln(iT)a 7 ^““/" 


'K 


e maXjj<|'£CT^rs'] Pn < Fi] 


^ / —«/3 

< exp ( 
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where the last inequality holds, because we can apply Proposition |9^ 


max 

n< [ecrx-^l 




'^n\j\erTKK^ ^olTo>r|-j^^K-| ] 


n<\eaKK^ P^^Tq > T\eaKK^\ 


(9.41) 


< 0(ln(K)cT 




On the other hand, we have 


^\eaKK^ [T\eK^ > Tq] = 1 - ^ ^ 


(d/ 6 )h^l - 1 


since d/h = 1 — 0{aK) and ax^K. 


(9.42) 

□ 


Proposition 9.5. Let (Zf) n>o a sequence of discrete time Markov Chain with state space Z 
and with transition probabilities 


\ + Cax, ifj = i + l, 


(9.43) 


^[Zn+i =j\Zn=i]= Pi^J) = ^ i - Cax, ifj = i-l, 

0 , else, 

for some constant C 7 ^ 0. Let Ti be the first hitting time of level i by and let Fj denote the 
law of conditioned on Zq = j and let ax a zero sequence such that -c 1 . 

(a) If Z^ is slightly supercritical, i.e. C > 0, then, for all i>l 

lim exp(ii'") Pi|-(^/2)f7xiri [^i-l)\{e/ 2 )(JKK^ <'r(i+i)r(e/2)(7xiri] =0- (9.44) 


K^oo 


(b) If Z^ is slightly subcritical, i.e. O' < 0, then, for all constants Ci, C 2 , Cs > 0 

lim exp(ii'") P(Ci+C 2 )re(TK^s:l [’'{Ci^C 2 +Ci)\LaKK^ <TCi\eaKK-\] = 0. 


K^oo 


(9.45) 


Proof. Since the transition probabilities of Z^ do not depend on the state of Z^, we have that 

^i\{e/2)<TKK^ [T{i-l)\{e/2)aKK^ > ^i+l)\{e/2)aKK^\ = ^\{e/2)aKK^ [ro > T2\{e/2)aKK^] (9.46) 


By ( |9.14| ) the left site of ( |9.46[ ) is equal 

1 - (1 - 2Cax + 0(a|))r("/2)-7Kiri 


1 - (1 - 2C'ax + 0(a|))2r(^/2)-xi^l 


> 1 — exp(—iP^a), 


since <JxK 3> With the same arguments, we obtain also (9.451. 


(9.47) 

□ 
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